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ABSTRACT 


Nonlinear six degree-of-freedom simulations for a variety of aircraft were created 
using MATLAB. Data for aircraft geometry, aerodynamic characteristics, mass / inertia 
properties, and engine characteristics were obtained from open literature publications 
documenting wind tunnel experiments and flight tests. Each nonlinear simulation was 
implemented within a common framework in MATLAB, and includes an interface with 
another commercially-available program to read pilot inputs and produce a three- 
dimensional (3-D) display of the simulated airplane motion. Aircraft simulations include 
the General Dynamics F-16 Fighting Falcon, Convair F-106B Delta Dart, Grumman F-14 
Tomcat, McDonnell Douglas F-4 Phantom, NASA Langley Free-Flying Aircraft for Sub- 
scale Experimental Research (FASER), NASA FFL-20 Lifting Body, NASA / DARPA X- 
31 Enhanced Fighter Maneuverability Demonstrator, and the Vought A-7 Corsair II. All 
nonlinear simulations and 3-D displays run in real time in response to pilot inputs, using 
contemporary desktop personal computer hardware. The simulations can also be run in 
batch mode. Each nonlinear simulation includes the full nonlinear dynamics of the bare 
airframe, with a scaled direct connection from pilot inputs to control surface deflections 
to provide adequate pilot control. Since all the nonlinear simulations are implemented 
entirely in MATLAB, user-defined control laws can be added in a straightforward 
fashion, and the simulations are portable across various computing platforms. Routines 
for trim, linearization, and numerical integration are included. The general nonlinear 
simulation framework and the specifics for each particular aircraft are documented. 



NOMENCLATURE 


a angle of attack, deg or rad 

ft angle of sideslip, deg or rad 

At numerical integration time step, s 

8 ' a , i e ft left aileron control surface deflection, deg 

8 a right right aileron control surface deflection, deg 

8 a aileron control surface deflection = [S a right - 8 a ie f t )^2 , deg 

8 C canard control surface deflection, deg 

S df differential flap deflection, deg 

S e elevator control surface deflection, deg 

Sy leading edge flap deflection, deg 

S lg landing gear deflection, deg 

S nf negative flap deflection, deg 

S pf positive flap deflection, deg 

8 r rudder control surface deflection, deg 

8 s b speed brake deflection, deg 

8 sp symmetric wing spoiler deflection, deg 

8y trailing edge flap deflection, deg 

<>th throttle deflection 

8 th i e f t left engine throttle deflection 

8 th, right right engine throttle deflection 

A outboard thrust-vectoring paddle deflection, deg 

0 Euler roll angle, rad 
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r 

0 


flight path angle, rad 
Euler pitch angle, rad 
p atmospheric density, slug/ft 3 

T eng engine power level time constant, s 

¥ Euler yaw angle, rad 

(O eng angular velocity of the engine rotating mass about the X body-axis, rad/s 

a x linear acceleration along X body-axis, g 

a Y linear acceleration along Y body-axis, g 

a z linear acceleration along Z body-axis, g 

A system matrix 

b wing span, ft 

B input matrix 

c mean aerodynamic chord, ft 

C output matrix 

c vector of mass / inertia properties 

C1...C9 inertia constants 

Ca non-dimensional axial force coefficient, C A = -C x 

Cd drag coefficient 

Cn minimum drag coefficient 

min 

Ci non-dimensional rolling moment coefficient 

Cl lift coefficient 

C j lift-curve slope 

C m non-dimensional pitching moment coefficient 

Cn non-dimensional normal force coefficient, C N = -C z 

C n non-dimensional yawing moment coefficient 

C x non-dimensional X body-axis force coefficient 

Cy non-dimensional Y body-axis force coefficient 
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Cz non-dimensional Z body-axis force coefficient 

D feed-through matrix 

2 

g earth gravitational acceleration = 32. 174 ft/s 

h altitude above mean sea level, ft 

Kng engine angular momentum about X body-axis, slug-ft'/s 

2 

I eng X body-axis moment of inertia for the engine rotating mass, slug-ft" 

2 

I x X body-axis moment of inertia, slug-ft 

I xz X-Z body-axis product of inertia, slug-ft 

I Y Y body-axis moment of inertia, slug-ft 2 

I z Z body-axis moment of inertia, slug-ft" 

m aircraft mass, slugs 

M Mach number 

O ( ) order 

p X body-axis angular velocity component, rad/s 

P c commanded engine power level, percent 

P a actual engine power level, percent 

q Y body-axis angular velocity component, rad/s 

q dynamic pressure, lb/ft 

r Z body-axis angular velocity component, rad/s 

2 

S wing reference area, ft 

T engine thrust, lb 

T temperature factor 

Tme idle engine thrust, lb 

T/ rolling moment about X body-axis generated by thrust, ft-lb 

T m pitching moment about Y body-axis generated by thrust, ft-lb 

T max maximum engine thrust, lb 

T m a X , a b engine thrust with maximum afterburner, lb 
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T mi a military engine thrust, lb 

T m in,ab engine thrust with minimum afterburner, lb 

T n yawing moment about Z body-axis generated by thrust, ft-lb 

T x X body-axis thrust component, lb 

T y Y body-axis thrust component, lb 

T z Z body-axis thrust component, lb 

u X body-axis translational velocity component, ft/s 

u control vector 

v Y body-axis translational velocity component, ft/s 

V t airspeed, ft/s 

w Z body-axis translational velocity component, ft/s 

x state vector 

x c .g. longitudinal center of gravity location, fraction of c 

x c .g.,ref reference longitudinal center of gravity location, fraction of c 
x E X earth- axis position, ft 

y output vector 

y E Y earth-axis position, ft 


SUPERSCRIPTS: 

time derivative 

SUBSCRIPTS: 

0 basic or trim value 

left left 

right right 

S stability axes 

E earth axes 


x 



NON-DIMENSIONAL DERIVATIVES: 


Longitudinal 
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dC z 
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c L 

q d qc 

2V t 
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2V t 
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Lateral /Directional 
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I. INTRODUCTION 


Nonlinear aircraft simulations are useful for dynamic analysis, control law design and 
validation, guidance and trajectory studies, air combat investigations, pilot training, and 
many other tasks. Most real-time nonlinear aircraft simulations are implemented in 
programming languages such as FORTRAN or C++, for historical reasons or to achieve 
the required execution speed. Advances in desktop computing capabilities have made it 
possible to implement real-time nonlinear aircraft simulations in an advanced computing 
environment called MATLAB 1 , using relatively inexpensive desktop personal computer 
hardware. MATLAB is easy to leam and program, is portable across many computing 
platforms without modification, and is a popular tool for scientific and engineering 
computations. The speed of current desktop computing hardware is sufficient for real- 
time execution of nonlinear aircraft simulations programmed in MATLAB, including 
numerical integration of the nonlinear equations of motion, general nonlinear 
aerodynamic and engine model computations, and interface with software to read pilot 
inputs and display of a pilot view or a 3-D view of the aircraft in motion from an 
arbitrary vantage point outside the aircraft. 

In addition to these computational hardware and software capabilities, reports 
available in the open literature contain aircraft geometry, mass / inertia properties, and 
extensive wind tunnel and flight test aerodynamic data for a variety of different aircraft, 
many of which are still in service. Data on engine perfonnance is less common in the 
open literature, but available nonetheless. 

This report describes realistic and easy-to-use nonlinear aircraft simulations 
implemented in MATLAB for a variety of different aircraft whose particular 
characteristics were available in the open literature. The nonlinear simulations include a 
general framework for the nonlinear equations of motion and a set of utility programs 
that apply to all the aircraft. For each individual aircraft, routines that implement the 
specific aerodynamic, propulsion, and mass / inertia characteristics of the aircraft were 
developed and documented. 

Aircraft simulations included in this report are the General Dynamics F-16 Fighting 
Falcon, Convair F-106B Delta Dart, Grumman F-14 Tomcat, McDonnell Douglas F-4 
Phantom, NASA Langley Free-Flying Aircraft for Sub-scale Experimental Research 
(FASER), NASA HL-20 Lifting Body, NASA / DARPA X-31 Enhanced Fighter 
Maneuverability Demonstrator, and the Vought A-7 Corsair II. Complete MATLAB 
simulations for each individual aircraft were developed in independent software 
packages. Each software package includes the capability to calculate aircraft trim 
conditions, linearize the equations of motion, run batch simulations of linear and 
nonlinear aircraft dynamics, and perfonn real-time interactive flight simulation. A 
commercially-available program called Aviator Visual Design Simulator (AVDS) 2 is 
used to interface each MATLAB nonlinear aircraft simulation with a human pilot in real 
time. Software to run each MATLAB simulation and exchange data with AVDS in real 
time is included in the simulation packages, along with 3-D model definition files for 


1 



each aircraft. The 3-D model definition files are used by AVDS to create the 3-D display 
of the aircraft in flight. 

Full documentation of each nonlinear aircraft simulation is included in this work, 
along with instructions for implementing the simulation software. Examples illustrate the 
steps required to execute the trim, linearization, and nonlinear simulation routines, and 
run the MATLAB nonlinear aircraft simulations in real time using AVDS. A short 
tutorial detailing the creation and importation of a 3-D Computer-Aided Design (CAD) 
aircraft model into AVDS is also provided. 


II. EQUATIONS OF MOTION 

Aircraft nonlinear equations of motion were implemented in MATLAB. Nonlinear 
aerodynamic and propulsion models were used to simulate the applied forces and 
moments for each aircraft. Usually, the aerodynamic and propulsion models were 
implemented as linear interpolations using look-up tables in a database. In some cases, 
analytic polynomial models were used to model the aerodynamic and propulsion 
dependencies on state and control variables. Some simplifying assumptions were made 
in the development of the nonlinear aircraft simulations, as described next. 


ASSUMPTIONS 

Each modeled aircraft was assumed rigid with constant mass density and symmetry 
about the X-Z plane in body axes. Forces and moments acting on the aircraft were from 
aerodynamics, propulsion, and gravity. Thrust was assumed to act along the X body-axis 
and through the center of gravity. A stationary atmosphere was assumed, with aircraft 
flight limited to altitudes less than 50,000 feet and subsonic speeds. The earth’s 
curvature was ignored and the earth was assumed to be fixed in inertial space, so that 
earth axes were considered inertial axes. The gravity field was assumed to be uniform, so 
that the center of mass and the center of gravity were coincident, and there were no 
gravity moments or changes in gravity with altitude. 

Full nonlinear six-degree-of-freedom rigid body translational and rotational aircraft 
motion was modeled, neglecting relative motion of the aircraft internal components, 
structural distortion, and sloshing of liquid fuel. The gyroscopic effects caused by 
rotating engine turbomachinery were included for some aircraft. Airplane rotation 
conventions follow the usual right-hand rule. Control surface deflections follow a 
standard right-hand rule sign convention, in which the fingers curl in the direction of 
positive control surface deflections when the thumb is pointed along the control surface 
hinge line in the direction of a positive body axis. 


COLLECTED EQUATIONS OF MOTION 

The six nonlinear rigid-body equations of motion were used to model the nonlinear 
aircraft dynamics in translational and rotational motion. These equations are 3-7 : 


2 



u =rv-qw- g sin 0 + 


qSC x +T 


m 


v = pw - ru + g cos 0 sin & + 


qSC Y 

m 


w = qu- pv + g cos 0 cos & + 


qsc z 

m 


(i) 


pi X~ r I xz = qSbCi -qr{l z -Iy) + qpl xz 

qI Y =q ScC m ~ pr (I x -I z )-(p 2 -r 2 )l X Z~ rh eng ( 2 ) 

rI Z ~ M XZ =qSbC n - pq{l Y -Ix)~ c l rI XZ + VKng 


where h eng is the magnitude of the angular momentum vector for the rotating mass of 
the engine, assumed to act along the positive X body-axis, 


heng 


heng 0 0 


I eng® eng ^ ^ 


.1? 


It is preferable to write the translational equations in terms of V t , cc, and jB instead of 
u,v, andw, because V t ,a, and jB can be measured directly on real aircraft, and have a 
more direct relationship to piloting and the aerodynamic forces and moments. The 
relationships between V t ,cc,(3 and u,v,w are: 


and 


V t = sfu 2 + v“ + w 


-1 


a = tan 


P = sin 1 


2 -f ' " 2 

(-] 

K u) 

( \ 

V 

Jt) 


u = V t cos a cos p 
v = V t sin P 
w = V t sin a cos P 


( 3 ) 


(4) 


Differentiating the equations (3) with respect to time gives 
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V t = 


a = 

fi = 


uu + vv + ww 

V t 

uw-wu 

2 . 2 
u + w 


V t v~vV ( 



1 

( \ 

2 

jr2 

r 'i 

h- 

V 



U ) 



( 5 ) 


The translational states for the nonlinear aircraft simulations are V t , (X, and (3 . The 
time derivatives of these translational states are found by computing u, v, and w from (1), 
then computing u, v, and w from (4) using current values of V t , a, and f3 , then 
substituting these u, v, and w values and the current value of V t into (5) to obtain 
V t , a, and j3 for use in the numerical integration routines. This approach avoids the 

lengthy nonlinear calculations associated with the equations of motion written in wind 
axes, and allows aircraft non-dimensional aerodynamic coefficient data to be 
implemented in the body-axis coordinate system, with no need for converting the data to 
wind axes components. 

For the rotational motion, straightforward algebraic manipulation transforms the 
equations (2) into a form suitable for numerical integration, with only one time derivative 
on the left side of each equation, 

P = (c\r + c 2 p + c 4 h eng )q + qSb (c 3 C 7 + c 4 C„ ) 

q = [c 5 p - c 7 h eng ) r - c 6 [p 2 - r 2 ) + qS cc 1 C m (6) 

r = (c 8 P - c 2 r + c 9 h eng )q + qSb (c 4 C l + c 9 C n ) 
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where 


c - 1 = 


{h ~b)b ~ bz' 
bb ~ 1 xz 1 


Cl = 


{Jx_ ~ l Y + l Z ) 1 XZ 

bb - bz~ 


Ci =■ 


bb bz 


c 4 


'AZ 


^x^z bz 


C 5 


^z ~b 


c 6 




(V) 


c 7 - ■ 


C 8 = 


( b ~b)b~bz ~ 


c 9 = 


'x 


bb bz 


The nonlinear aircraft simulation also includes rotational kinematic equations and 
navigation equations. The rotational kinematic equations, which relate Euler angular 
rates to body-axis angular rates, are given by: 

<P = p + tan0[qsin0 + rcos0 ) 

0 = q cos 0 — r sin O (8) 

K p _ q sin <P + r cos 0 
cos 0 


Equations (8) are nonlinear state equations for the Euler angles 0, 0, and W , in a 
form suitable for numerical integration. The rotational kinematic equations describe the 
time evolution of the aircraft attitude angles, which are required to properly resolve the 
gravity force along the aircraft body axes in (1). A singularity exists for the rotational 


kinematic state equations at 0 = ±90° . This is eliminated within the simulation code by 


limiting 0 to the ranges 0° to ±89.99° and ±90.01° to ±180°. This limits 


(cos@) 1 


to a 


maximum value of approximately 5730, which allows sufficiently accurate calculations 
of the kinematic state derivatives near the singularity. In the nonlinear simulation, 0 
automatically skips from ±89.99° to ±90.01° and vice-versa depending on the direction of 

aircraft pitch angle motion near 0 = ±90° . 


The navigation equations relate aircraft translational velocity components in body 
axes to earth-axis components, neglecting wind effects. These differential equations 
describe the time evolution of the position of the aircraft c.g. relative to earth axes: 
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( 9 ) 


Xfr = u cos F cos 0 + v (cos F sin 0 sin 0 - sin F cos 0) 

+ w ( cos F sin 0 cos 0 + sin F sin0) 
y E = u sin F cos 0 + v(s in F sin 0 sin 0 + cos F cos 0) 

+w(sin F sin 0 cos 0 - cos F sin0 ) 
h = u sin 0 - v cos 0 sin 0 - wcos 0 cos 0 

Assuming thrust acts along the X body-axis, body-axis accelerations a x , ciy ’ anc ' a z 
are calculated from 


_qSC x +T 

a x ~ 

mg 

qSC Y 

a Y = — 

mg 


qsc z 

mg 


( 10 ) 


Equations (1), (4), (5), (6), (7), (8), and (9) are the nonlinear aircraft equations of 
motion for each aircraft simulation. These nonlinear dynamic equations are implemented 
for each simulated aircraft in functions called ACname deq.m, where ACname is the 
aircraft name (e.g., F16_deq.m, F106B_deq.m, etc.). Body-axis accelerations 
a x , a Y , and a z are also calculated in each ACname_deq.m function, using equations 
(10). 

ENGINE MODEL 

The engine power dynamic response was modeled with an additional state equation 
as a simple first order lag in the actual power level response to commanded power level: 

Pa=—(P C ~Pa) (ID 

^ eng 

Commanded power level was computed as a function of throttle position, 

Pc = PAX) 
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and the engine power level time constant T eng was chosen at or near 1 sec. Engine thrust 

force was typically computed by linear interpolation of an engine thrust database, as a 
function of actual power level, altitude, and Mach number, 

T = T(P a ,h,M ) 

For each aircraft simulation, calculation of the engine thrust occurs in a function 
called ACname_engine.m, where ACname is the aircraft name (e.g., F16 engine. m, 
F106B_engine.m, etc.). The routine tgear.m implements the throttle gearing, which 
translates the throttle deflection in the interval [0,1] to commanded power level in the 
interval [0,100], but not always in a linear fashion. The routine rtau.m computes the 
engine power level time constant r eng , and pdot.m implements the engine thrust 

dynamics in equation (11). 

For many of the simulated aircraft, engine data was not available. In these cases, the 
engine model for the F-16 aircraft was used, and the output thrust was scaled to match the 
engine thrust of the specific aircraft. 


AERODYNAMIC MODEL 

Aerodynamic forces and moments acting on each aircraft were characterized in tenns 
of non-dimensional coefficients in body axes. Typically, the non-dimensional force and 
moment coefficients were calculated from an expression that included separate build-up 
functions for each non-dimensional coefficient. The value for each build-up function was 
usually found by linear interpolation of aerodynamic data, using current values of the 
states and controls. For values of states and controls outside the range of the available 
data, the interpolation routines extrapolated linearly using the nearest data points. In a 
few cases, build-up functions or complete non-dimensional coefficients were computed 
using analytical polynomial expressions involving state and control variables. Once the 
non-dimensional aerodynamic force and moment coefficients were calculated, these 
values were used in equations (1) and (6) with thrust, mass / inertia properties, geometry, 
and current states and controls to compute state time derivatives for the aircraft 
dynamics. Subsequent sections describing the individual aircraft simulations contain 
details of the aerodynamic model for each nonlinear aircraft simulation. 


ATMOSPHERIC PROPERTIES 

Air density and the speed of sound were calculated using relations modeling the U.S. 
Standard Atmosphere 3 . Required quantities that depend on these atmospheric properties, 
namely Mach number M, and dynamic pressure q, were also calculated. The 
relationships are: 


T = 1 - 0 . 703 x 10“ 5 /2 
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/ *\ 4.14 

p = 0. 002377 (r ) 


V, 


J 1.4(1716.3)390 


M = < 


V t 


L.4(l716.3)(519r*) 


h > 35,000 ft 
h< 35,000 ft 



( 12 ) 


The function atm.m was used to compute all the quantities from equations (12), given 
h and V t . 


AIRCRAFT STATES AND CONTROLS 

For a typical aircraft, the state vector was: 

x = [V t a (3 & & T p q r x E y E h P a ] T (13) 

The elements of the state vector are associated with equations (1), (4), (5), (6), (7), (8), 
(9), and (1 1). A typical control vector was: 

U = [ S lh S e S a S r] T ( 14 ) 

The number of control inputs and engine power states varied depending on the number of 
control surfaces and aircraft engines. The set of coupled, nonlinear, first-order ordinary 
differential equations that comprise the simulation model can be represented by the 
vector differential equation: 


x = f(x,u) (15) 

The output equations can be represented by the vector equation: 

y = h(x,u) 
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III. GENERAL SIMULATION TOOLS 


Approximating the solution to the nonlinear equations of motion requires a suitable 
numerical integration technique. Linearization of the nonlinear equations of motion 
requires appropriate numerical gradient techniques, which are necessary to calculate 
elements of the linearized system matrices. Numerical algorithms for both these tasks are 
included in the software. Reference [8] provides background infonnation on the 
numerical techniques applied. Runge-Kutta and Adams-Bashforth numerical integration 
routines are used to approximate solutions to the nonlinear equations of motion. Finite 
difference numerical gradient procedures are used to determine the elements of the 
linearized system matrices. The scripts called gen_ACname_m.odel.rn, where ACname is 
the simulated aircraft name (e.g., gen_F16_model.m, gen_F106B_model.m, etc.), 
demonstrate the use of the trim and linearization tools for each simulated aircraft. 


TRIM 

A linearized form of the state-space model for the nonlinear equations of motion is 
obtained using a first-order multivariable Taylor series expansion about an equilibrium 
flight condition known as trim. At the trim condition, all translational and rotational 
accelerations are equal to zero. This condition is described by the equation: 


x Q =0 = f(x 0 ,u 0 ) 


(16) 


where only the state equations for V t , a, (3, p, q, and r are included in equation (16). The 

six force and moment equations are augmented with a rate of climb constraint, and a 
coordinated turn constraint, for a total of eight trim equations 3 ’ 4 : 


0 = V t = 


0 = a = 
0 = J3 = 


uu + vv + ww 

V t 

uw-wu 

2 . 2 
u + w 


VtV-vVt 



1 

( \ 

2 

’1 

1 1 - 

\ V 



U j 
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( 17 ) 


0 = P = (c\r + c 2 p + c 4 h eng )q + qSb (c 3 C) + c 4 C „ ) 
0 = q = ( c 5 p - c 7 h eng ) r - c 6 [p 2 - r 2 ) + qS cc 1 C m 
0 = r = (c 8 p - c 2 r + c 9 h eng ) q + qSb(c 4 C l + c 9 C n ) 


^ ^ ab + sin yJ a 2 - sin y + b 2 

0 = tan 0 - 

a" - sin y 

0 = Geos f5 [sin a tan 0 + cos a cos O) 
where 


0 ^ 90° [rate of climb constraint] 

0 -t 90° [coordinated turn constraint] 


a = cos a cos [3 

b = sin O sin jB + cos 0 sin a cos [3 



g 


With eight trim equations to find the trim condition, any combination of eight state or 
control surface deflection variables can be free to vary so that the trim equations can be 
satisfied. 

User-defined values for fixed states, control surface deflections, and trim values for y 
and turn rate W are set in the function ACname trm.m, where ACname is the simulated 
aircraft name (e.g., F16_trm.m, F106B_trm.m, etc.), via code modification. The 
ACname trm.m function calls the ACname_deq.m function to implement the state 
derivative equations included in the trim equations. Pitch angle 0 is limited to the ranges 
0° to ±89.99° and ±90.01° to ±180° in the function ACname_deq.m, as described above. 

This eliminates the singularity at 0 = ±90° within the rate of climb and coordinated turn 
constraint equations. In addition, the engine power level dynamics in equation (11) are 
bypassed for the trimming operation, so P a is set equal to P c inside ACname trm.m. 

Trim solutions are obtained by solving the eight nonlinear trim equations above, 
using the function solve. m. Function solve. m uses a Newton-Raphson iterative solution 
technique, where gradients are calculated in grad.m by central finite differences. Trim 
solution convergence is determined using cnvrg.m. Convergence is achieved when the 
absolute value of every trim equation is less than the solution tolerance, and the absolute 
value of every change in the free states and controls for successive iterations has reduced 

_ o 

to less than the solution tolerance. The solution tolerance is specified as 1x10 within 
the code. When the trim search is successful, the words CONVERGENCE CRITERIA 
SATISFIED are output to the MATLAB command window. 
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LINEARIZATION 

A linear, time-invariant (LTI) system representing aircraft dynamics for small 
perturbations about a reference trim condition is given by the state and output equations: 


x = Ax + Bu 
y = Cx + Du 


(18) 


Elements of the constant-coefficient matrices in the linearized dynamic equations are 
gradients evaluated at the trim values x 0 and u {] : 




C = 


df(x,u) 


dx 

dh( x, u ) 
dx 


l 0 '“0 


n x n matrix 


/ x n matrix 


B = 


D = 


df(x,u) 


du 

dh( x, u ) 
du 


l o “o 


x Q ,u Q 


■ n x m matrix 


: / x m matrix 


(19) 


where n is the number of elements in the state vector, m is the number of elements in the 
control vector, and / is the number of elements in the output vector. The gradients are 
obtained numerically by perturbing each state and control input independently and 
recording the changes in the trimmed state and output equations. This is done using the 
numerical technique of central finite differences: 


df(x,u) 
dx 


1*0 ,H 0 


Ax , Ax 

f( x o+^’ u o)-f( x o-^’ u o) 

Ax 


df(x,u) 
du 


'*o.«o 


. , Au . . . Au , 

f( x o^ 0 +—)-f(x 0 ,u 0 - —) 

Au 


( 20 ) 


where Ax and Au are the perturbations of the state and control variables, respectively. In 
the simulations, the output vector is the state vector [h(x,u) = x) , so C is an identity 

matrix, and D is a null matrix. A flowchart describing the numerical linearization 
scheme is shown in figure 1 . 
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Figure 1. Flowchart for Numerical Linearization 

For a reference trim condition of straight and level flight, linearization results in two 
decoupled sets of linear, constant-coefficient differential equations for longitudinal and 
lateral / directional motion. The linearization is valid for small perturbations about the 
reference trim condition. 

Detailed Linearization Example 

The trim state and control vectors (xo and «o) are assembled in the function icjtrm.m, 
for use in Inze.m, which computes the linear system matrices A, B, C, and D, using 
central finite differences. Perturbations equal to 1% of the trim reference values are used 
for each state and control variable during gradient calculation for the elements of the 
system matrices. If the trim state or control is zero, the perturbation size is set to 0.01. 
Setting the variable lonflg to 1 within the MATLAB workspace before executing 
gen_ACname_model.m, where ACname is the simulated aircraft name (e.g., 
gen_F16_model.m, gen_F106B_mode/.m, etc.), produces the longitudinal system 
matrices. Setting lonflg to 0 provides the lateral / directional system matrices. 
Typical constant-coefficient longitudinal or lateral / directional system matrices A, B, C, 
and D, for an aircraft with controls for throttle, elevator, aileron, and rudder deflection, 
are created in the form: 
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For lonf lg = 1: 




~Wt~ 




~^v t ~ 

Aa 


Aa 


A8 t h 


Aa 

Aq 

= A 

Aq 

+ B 

A8 e 

y = c 

Aq 

A© 


A0 




A0 

1 

h, 

1 


1 

1 




1 

X] 

1 


+ D 


A8 th 

AS e 


For lonf lg ^ 1: 


~A(5~ 


~Ap- 




~AP 

Ap 

= A 

Ap 

+ B 

1 

1 

o 

II 

Ap 

Ar 


Ar 


L4J 


Ar 

A0 


AO 




AO 


A8 r 


The following example illustrates the procedure used to trim the F-16 simulation and 
calculate the longitudinal and lateral / directional linearized system matrices. 

Fixed Quantities Free States / Controls 

V t = 500 ft/s x E = 0 ft 

p = 0 rad/s y E = 0 ft 

q = 0 rad/s h = 10, 000 ft 


a, j8, 0, 0, 8ti„ 8 e , 8 a , and 8, 

r = 0 rad/s 7 = 0.349 

¥ = 0 rad ¥ = 0.052 rad/s 


x c g =0.3 

Step 1. Modify F16_trm.m to account for the desired trim condition, indicating the fixed 
and free state and control variables. Within this file, the variables xfree and 
xinit control how the states are handled for trim. The definitions of xfree 
and xinit are located on lines 117 and 118, respectively. Line 117 and 118 
must be modified to read: 

xfree = [0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0] ^ ; 
xinit = [500,0,0,0,0,0,0,0,0,0,0,10000,0]'; 

In the vector xfree, a 0 element indicates that the corresponding state is 
fixed, while 1 specifies that the corresponding state is free to be chosen by the 
trim solver. Values for the fixed states are input using the vector xinit. If the 
state is free, the corresponding element of xinit contains the starting value 
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used to determine trim solution. If the state is fixed, the corresponding element 
of xinit contains the fixed value. 

The variables ufree and uinit contain information for fixed and free 
control variables. They are located on lines 122 and 123, respectively. All 
controls are free variables for this particular example. Lines 122 and 123 should 
be changed to: 

ufree = [ 1 , 1 , 1 , 1 ] ; 

uinit = [ 0 , 0 , 0 , 0 ] ; 

Other starting values for the free state and control variables may need to be 
tried in order to converge to the trim solution. 

Flight path angle yis input in radians on line 127: 

gam=0 .349; 

which corresponds to an angle of 20 deg. 

Turn rate is input in rad/s on line 131: 

trate =0.052; 

This corresponds to a turn rate of 3 deg/s. This completes the necessary 
modifications to F16_trm.m. Save the modified file. 

Step 2. Modify the file F16_massprop.m to account for an x c . g . location at 30% of the 
mean aerodynamic chord. Percentages must be normalized to values between 0 
and 1 . Change line 45 of the file to read: 

xcg=0 . 3 ; 

Step 3. Set the variable lonflg in the MATLAB workspace to either 1 or 0. Setting 
lonflg to 1 will provide the longitudinal system matrices, while setting the 
variable to 0 will provide lateral / directional matrices. 

Step 4. Run the file gen_Fl 6_model.m by typing the following on the MATLAB 
command line: 

gen_Fl 6_model ; 

This function computes the trim solution as well as the linear system matrices at 
the desired flight condition and longitudinal center of gravity location. After 
program execution, x 0 and u 0 are stored in the MATLAB workspace as 

variables xO and uO, respectively. The linear system matrices are stored in the 
MATLAB workspace as matrix variables A, B, C, and D. 


The preceding procedure can be applied to the other supplied aircraft simulations with 
little modification (e.g., replace F16 with F106B, F14, F4 , etc. in the called functions). 
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NUMERICAL INTEGRATION 


The equations of motion for a rigid aircraft are a set of coupled first-order nonlinear 
ordinary differential equations. Experimental aerodynamic and thrust data are used to 
model the applied aerodynamic and propulsion forces and moments for arbitrary states 
and controls. There is no closed form solution to such problems, so the equations must 
be solved using numerical integration. Techniques for solving this initial value problem 
for ordinary differential equations are employed to obtain approximate solutions at 
discrete points along the aircraft state trajectory. 

The equations of motion have already been arranged into the fonn: 

x = f(x,u) (21) 


where 


x = 


V t a fi 0 0 W p q r x E y h P a 


« = [ 4 , 4 4 4 Y 


( 22 ) 


Control vector elements can be added or removed depending on the modeled aircraft. 

A number of numerical integration techniques are available to approximate the 
solutions of the well-posed equations of motion. The choice of which technique to use 
depends on the desired accuracy of the solution and the computational effort to be 
expended. Three types of numerical integration techniques are employed within the 
software to approximate solutions to the equations of motion - Euler’s method, 
Runge-Kutta method of orders two and four, and Adams-Bashforth method of order 
three. A flowchart illustrating the procedure used during the nonlinear simulation is 
given in figure 2. 
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Figure 2. Flowchart for Nonlinear Simulation 


Euler’s Method 

Euler’s method is a one-step integration routine that approximates solutions to the 
initial value problem at discrete time steps At. Solutions to the state equations are 
approximated starting from initial states given at t = 0 with arbitrary control inputs 
provided for t >0. The local truncation error for this method is O(At). A single state 
derivative calculation is required for each time step, making this the least 
computationally intensive of the three approximation methods used in the simulation. 
However, it is also the least accurate. The algorithm is given in difference equation form 
as: 


x(t + At) = x(t) + At f[x(t),u(t )) 


(23) 


Runge-Kutta Method of Order Two 

The second-order Runge-Kutta method is also a one-step integration routine but 
provides for a more accurate solution approximation as compared to Euler’s method. 
However, two state derivative calculations per time step are required making it 
comparatively slower than Euler’s method. The local truncation error for this routine is 
0(At~). The difference equation is given by: 


x(t + At) = x(t) + Atf 


\L 


x(t) + Yf{x(t)’U(t)) 



( 

,u 

t H 

l 2 JJ 


(24) 
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Runge-Kutta Method of Order Four 

The fourth-order Runge-Kutta method is a one-step routine with a local truncation 
error of 0(At 4 ). Solution approximations are more accurate than either the Euler or 
second-order Runge-Kutta method, but the method runs one-fourth and one -half as fast, 
respectively. The difference equation for this scheme is given by: 


x(t + At) = x(t ) + At 


^ k") k-, k/L 

— H — — + — -t — — 


v 


(25) 


where 


h = f{x(t),u(t)) 


f s At , 

' At\] 

x(t)-\ L,u 

t H 

7 

V z 

v 2 

f , At , 

( At\] 

x( t)- 1 k 2 ,u 

t H 

V 2 

V 2 


k 4 = f[x(t) + Atk 3 ,u(t + At)) 


Adams-Bashforth Method of Order Three 

The third-order Adams-Bashforth method is a multi-step routine that approximates 
the solution of the initial value problem using three previous step calculations. The 
approximation has a local truncation error of 0(At ). Second-order Runge-Kutta is used 
for the initial time steps before switching to the Adams-Bashforth integration method. 
The difference equation for the third-order Adams-Bashforth method is given by 


x(t + 3At) = x(t + 2 At) + ^^23 f [x(t + 2At),u(t + 3At)) 

-16 f[x(t + At),u(t + 2At)) + 5 f [x(t),u(t + At))~^ 


(26) 


Detailed Numerical Integration Example 

The state vector derivative x is integrated numerically in the function ACname.m, 
where ACname is the simulated aircraft name (e.g., F16.m, F106B.m, F4.m, etc.) using 
arbitrary control inputs for t > 0 with initial states specified at t = 0. This produces a set 
of state time histories x{t) and output time histories y(t) for / > 0. Inputs to ACname.m 
are the control input time history u(t), time vector t, initial state vector x at t = 0, and the 
vector of mass / inertia properties c. The function mksqw.m is provided to allow creation 
of a u(t) comprised of arbitrary square waves with user-specified amplitudes, widths, 
delay, step size, and total length in seconds. However, any u(t) can be used with the 
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nonlinear simulation. A user-specified integration routine, such as the Adams-Bashforth 
method of order three ( ab3.m ), Runge-Kutta method of order two ( rk2.m ), or Runge- 
Kutta method of order four ( rk4.m ), is called within the ACname.m function. 

A sample procedure for calculating state and output trajectories is provided below for 
the F-16 simulation. Procedures for the other aircraft simulations are similar. 

Step 1. Create u(t) along with t, which will be inputs to the nonlinear simulation 
function. The control vector must be stored in the MATLAB workspace as a 
matrix with the individual controls making up the columns and the deflection of 
each control over time comprising the rows. The time vector contains the 
discrete time points for the control vector. An example using the supplied 
function mksqw.m to create u(t) is shown below in figure 3. The elevator control 
surface 8 e will be deflected using a doublet with an amplitude of 1 degree that 
will be applied after a 3 second delay, while keeping the remaining control 
inputs {8th, 8 a , and 8,) unperturbed from their trim values. Each pulse of the 
doublet will have a width of 2 seconds. The entire length of the maneuver will 
be 10 seconds with a sampling interval 0.025 seconds. The elevator perturbation 
time history is shown in figure 3. 


0 
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O 

0 
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0 
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-1.5 
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Figure 3. Elevator Perturbation 


a. Type the following command in MATLAB: 

[usw, t ] =mksqw (1,2, [ 1 1], 3,0.025,10); 

The first input to the function gives pulse amplitude, the second gives 
the width of a single pulse in seconds, the third gives the vector of 
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integer pulse widths, the fourth gives the delay in seconds, the fifth 
gives the sampling time, and the final input gives the length of the entire 
maneuver in seconds. The output of the function is stored in the 
MATLAB workspace under the variable names usw and t. The 
variable usw is a column vector composed of 401 elements describing 
the deflection history. The variable t is also a column vector composed 
of 401 elements recording the sample times for usw. 

b. Create the control matrix u within the MATLAB workspace by adding 
usw to the second column of a 401 x 4 dimensional matrix, where the 
first column corresponds to values of 5th, the second to 5 e , the third to 5 a , 
and the last column to 5,-. If the elements in the columns of u are 
initially set to the trim values of the controls, then adding usw to 
column 2 creates u(t) for an elevator doublet perturbation deflection, 
leaving other control inputs unperturbed from trim values during the 
maneuver. 

Step 2. Define the initial state vector xO by creating a 13 x 1 vector composed of the 
user-defined state initial conditions. This vector contains the same state 
variables as the state vector x of the simulation model. To start from trim, use 
the trim state vector output from ic Jtrm.m or gen ACname model.m. 

Step 3. Generate the mass property vector c by running the function FI 6_massprop.m: 

c=F 1 6_mas sprop ; 

The vector c will be available in the MATLAB workspace after execution of this 
function. Center of gravity position x c . g , can be adjusted within this file to 
simulate various longitudinal c.g. locations. 

Step 4. Run the nonlinear F-16 simulation using the following command on the 
MATLAB command line: 

[y, x] =F16 (u, t, xO, c) ; 

The output time history is stored in the MATLAB vector y, and the state output 
time history is stored in the vector x. 


By default, the nonlinear simulation implements the Runge-Kutta integration method 
of order two. The user can modify the code to implement either the Runge-Kutta method 
of order four or the Adams-Bashforth method of order three. This is done by modifying 
the function F16.m following the procedure below: 

To use Fourth-Order Runge-Kutta 

Change line 102 to read: 

%x=rk2 ( ' FI 6_deq ' , u, t , xO , c) ; 

line 101 to read: 

x=rk4 ( ' FI 6_deq ' , u, t , xO , c) ; 
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and line 103 to read: 

%x=ab3 ( ' FI 6_deq ' , u, t , xO , c) ; 
To use Third-Order Adams -Bashforth 
Change line 102 to read: 

%x=rk2 ( ' FI 6_deq ' , u, t , xO , c) ; 
line 101 to read: 

%x=rk4 ( ' FI 6_deq ' , u, t , xO , c) ; 
and line 103 to read: 

x=ab3 ( ' FI 6_deq ',u,t,x0,c) ; 
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IV. AVDS / MATLAB INTERFACE 

A MATLAB software interface to the software package called AVDS was developed 
to provide the capability to fly the MATLAB nonlinear aircraft simulations interactively 
in real time, using AVDS to read pilot inputs and display the aircraft response. Pilot 
inputs are measured by AVDS and used as inputs to the MATLAB simulations. Results 
of the simulation numerical integration are sent to AVDS and used to draw 3-D images of 
the aircraft in motion. Control deflections are also shown on the 3-D aircraft model 
displayed by AVDS. The following material provides instructions for AVDS / 
MATLAB software installation and usage. 


EQUIPMENT / REQUIREMENTS 

The nonlinear simulation software described here was developed and tested on 
IBM-compatible personal computers (PC), using MATLAB version 6.1, and AVDS 
version 1.2.14. System requirements for these commercially-available software tools are 
specified in their documentation. A graphics accelerator card capable of hardware 
acceleration using OpenGL™ produces smoother-running piloted simulations. 

Control inputs are supplied to AVDS via a mouse, joystick, and / or keyboard. A 
joystick is recommended. A Microsoft® Sidewinder® Precision 2 joystick has been used 
exclusively during software development and is recommended for use; however, any 
joystick can be used. This particular joystick contains a throttle control lever, hat switch, 
rudder control, and eight programmable buttons, which provide for somewhat realistic 
aircraft controls. 


SOFTWARE INSTALLATION 

AVDS is available from RasSimTech Ltd. and is not part of the MATLAB aircraft 
simulation packages. The AVDS software allowing data communication between 
MATLAB and AVDS is available with version 1.2.14 of AVDS, available on the Internet 
at http://www.rassimtech.com . 

Once the AVDS software is installed, the 3-D aircraft images provided for each 
simulated aircraft can be installed into AVDS. The 3-D image files are included in the 
aircraft simulation packages. The provided 3-D aircraft models are installed into AVDS 
using the following procedure: 

Step 1. Copy the following aircraft 3-D files to the AVDS directory 

\usr\local\A VDS\craft\: 

\usr\local\F-l 6\F1 6_hires. txt 
\usr\local\F-16\F16_lores.txt 
\usr\loca/\F-106B \F1 06B. txt 
\usr\local\F-14\F14.txt 
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\usr\local\F-4\F4.txt 
\usr\locaI\FASER \FASER_hires. txt 
\usr\loca/\FASER \FASER_lores. txt 
\usr\local\HL-20\HL20.txt 
\usr\loca\XX-31XX31.txt 
\usr\locaIA--7\A 7. txt 

High and low resolution 3-D aircraft models are provided for the F-16 and 
FASER aircraft. The user has the choice to use either during interactive piloted 
simulations. Low-resolution models are provided for simulation on slower 
computers. 

Step 2. Add the following lines to the fde craftcap.txt located in the directory 
\usr\local\A VDS\craft\: 

include "F16_lores.txt" 

include "F16_hires.txt" 

include "F106B.txt" 

include "F14.txt" 

include "F4.txt" 

include "FASER_lores.txt" 

include "FASER_hires.txt" 

include "HL20.txt" 

include "X3 1 .txt" 

include "A7.txt" 

This file tells AVDS which 3-D aircraft models to make available for interactive 
simulation. 


INTERACTIVE SIMULATION 

The user can interactively provide control inputs through AVDS using a mouse, 
joystick, and/or keyboard. For control input device selection and configuration, consult 
section 3. 4.2. 9 of the AVDS User’s Manual". The following will discuss providing 
interactive control input with the Microsoft® Sidewinder R Precision 2 joystick. Other 
joysticks should work as well, but only this joystick was used in the development and 
testing of the nonlinear simulations. Configuration of the initial aircraft states is also 
described. 
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Joystick Control Input 

Control inputs for pitch, roll, yaw, and thrust during piloted interactive simulations 
can be provided by a Microsoft® Sidewinder* Precision 2 joystick. For proper set-up of 
the joystick, set the scale factors for X, Y, Z, and R joystick axes to 2, 2, 1, and -2, 
respectively, within the AVDS joystick configuration screen (refer to section 3. 4.2. 9 of 
the AVDS User’s Manual ). This nonnalizes AVDS inputs to the ranges 0 tol for thrust 
control, and -1 to 1 for pitch, roll, and yaw. The normalized values are then converted to 
actual control surface deflections, measured in degrees, within the function 
AVDS_Matlab_ACname.m, where ACname is the simulated aircraft name (e.g., 
AVDS_Matlab_F16.m, AVDS_Matlab_F106B.m, AVDS_Matlab_F4.m, etc.). Other 
joystick settings, such as dead-zone, bias, trim step size, etc., can also be configured 
within the AVDS joystick configuration screen. 

Aircraft pitch control is achieved with the fore / aft movement of the joystick, 
left / right movement of the joystick provides roll, and clockwise / counterclockwise twist 
of the joystick supplies yaw control. Throttle deflection is actuated by means of a small 
lever located at the base of the joystick. A hat switch located on the tip of the joystick is 
used to set elevator and aileron control surface trim deflections. Trim elevator deflection 
can be varied by pressing the top or bottom of the trim hat, while trim aileron can be set 
by pressing the left and right sides of the trim hat. The eight programmable buttons of 
the joystick have not been configured for use in AVDS / MATLAB interactive 
simulations. 

Initial States for Interactive Simulation 

AVDS allows specification of aircraft initial conditions for interactive simulations. 
The user can specify the aircraft starting altitude, velocity, Euler angles, and throttle 
position. These are set in the Aircraft State Initialization pop-up window, shown in 
figure 4, which is called within AVDS prior to simulation start-up. This is done as 
follows: 

Step 1. Select Initialize from the top menu 
Step 2. Select Aircraft State 


Initial altitude, given in feet above sea level, is adjusted by moving the corresponding 
slider up or down. Initial body-axis velocity components u, v, and w (in ft/s), initial Euler 
angles <P, 0, and l F (in degrees), and initial throttle position S t /, (in percent) are adjusted 
similarly. Starting xe and ys positions cannot be adjusted with these sliders. The default 
starting positions for xe and ve are both set to 0 ft. Before hitting the OK button to 
complete initial state configuration, make certain that the Apply Settings button is 
depressed so that a check appears in the box. This allows the initial states set with the 
sliders to take effect in the interactive simulation. The button Apply Init Lat/Long has 
no effect within the AVDS / MATLAB simulation since starting xe and ve aircraft 
positions are set within the MATLAB simulation software. 
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Figure 4. AVDS Aircraft State Initialization Window 

A VDS /MATLAB Piloted Simulation Software 

Piloted simulations are executed using the function AVDS_Matlab_ACname.m 
included in the aircraft simulation packages. This function establishes communication 
between AVDS and MATLAB and implements the integration routine to numerically 
solve the nonlinear equations of motion. Integration time steps during piloted 
simulations are dynamically updated within the AVDS_Matlab_ACname.m file, to 
prevent MATLAB from lagging behind AVDS simulation time. This allows 
synchronization of the real-time simulation of the aircraft in MATLAB with the 3-D 
display in AVDS. However, an upper limit (0.035 s) is imposed on the time steps to 
prevent loss of model fidelity. For less powerful computers (e.g., 500 MHz and slower), 
this can cause MATLAB to lag behind AVDS during simulation. 

Control inputs at a given time step are sent from AVDS to MATLAB via a shared 
memory connection. Solutions of the nonlinear equations of motion are computed in 
MATLAB using the control inputs from AVDS. The solution is returned to AVDS by 
MATLAB, using the same data pathway. Control inputs and aircraft initial states, 
supplied by AVDS, are stored in the MATLAB workspace using the variable name 
InputVector. The output vector, stored in the variable called OutputVector, 
contains the translational and rotational states, angle of attack, sideslip angle, linear 
accelerations in g’ s, airspeed, and Mach number, computed by the nonlinear simulation 
in MATLAB. Elements within InputVector and OutputVector are identified to 
AVDS via the file ACname.sim.ini, where ACname is the simulated aircraft name (e.g., 
F16.sim.ini, F106B.sim.ini, etc.). Refer to section 4.2 of the AVDS User’s Manual" for 
modification instructions for this file. Typical elements of these vectors are shown 
below: 
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InputVector = [flagrun 4 4 4 4 h <AVDS 0 0 h 0 0 ¥ V ( P ( , (X fif 
OuputVector = \jx, E y E hO0¥ a /3gV t M [control surface deflections]] 7 ^ 

where flag run is a flag equal to one when AVDS is running, and zero otherwise, and 
t A vDS ' s ^e simulation time in the AVDS program. 

Figure 5 shows a flowchart illustrating the interactive simulation. 



Figure 5. Interactive Simulation Flowchart 

The translational and rotational aircraft states calculated within MATLAB are used 
by AVDS to display the 3-D aircraft model along the calculated trajectory. Quantities 
important to the pilot ( a , /?, M, V t , h, 0, and a N (g’s)) are displayed to the user in the 

AVDS Head-Up Display (HUD). A screen shot of the AVDS screen during interactive 
simulation of the F-16 is shown in figure 6. Flight path of the aircraft is shown with 
trailing ribbons generated on each wing tip and the aircraft center tail section. 
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Figure 6. Screen Shot of the F-16 AVDS / MATLAB Interactive Simulation 


Stick Gain Modification 

Stick gains for elevator, aileron, and rudder are applied within the function 
AVDS_Matlab_ACname.m to allow for easier open-loop flying. All stick gains for the 
simulated aircraft are set to 0.5 by default. The procedure used to modify the stick gains 
for the F-16 simulation are given below (other simulation modifications are similar): 

To change elevator stick gain 

Modify line 123 within AVDS_Matlah_F16.m to read: 

gain_el = xxx; 

where xxx is the user-defined numerical value for the elevator stick gain. 

To change aileron stick gain 

Modify line 124 within AVDS_Matlab_F16.m to read: 

gain_ail = xxx; 

where xxx is the user-defined numerical value for the aileron stick gain. 

To change rudder stick gain 

Modify line 125 within AVDS_Matlab_F16.m to read: 

gain_rdr = xxx; 

where xxx is the user-defined numerical value for the rudder stick gain. 
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Control Deflection Exaggeration Factor 

Deflections of aircraft control surfaces, shown in green on the AVDS aircraft 3-D 
models, can be visually exaggerated for easier viewing. By default, all deflections shown 
in AVDS are exaggerated by a factor of two. The following provides an example of 
modifying the exaggeration factor for the F-16. Other simulation modifications are 
similar. 

To adjust control surface deflection exaggeration factor 

Modify line 130 within AVDS_Matlab_F16.m to read: 

cdfac = xxx; 

where xxx is the user-defined numerical value for the deflection exaggeration 
factor. If cdfac is set to 2, the displayed control surface deflections will be 2 
times the actual control surface deflections being used in the MATLAB 
nonlinear aircraft simulation. 

Integration Routine 

Integration routines used to solve the equations of motion for the AVDS / MATLAB 
interactive simulations are implemented within AVDS Matlab ACname.m as opposed to 
ab3.m, rk2.m, and rk4.m for the batch nonlinear simulations. The 
AVDS_Matlab_ACname.m function can implement either the Adams-Bashforth 
integration method of order three, the Runge-Kutta method of order two, or Euler’s 
method. The following gives an example for implementing particular integration 
routines for the F-16 simulation (. AVDS_Matlab_F16.m ). Modifications for other 
simulations are similar. 

To implement Adams-Bashforth method of order three 

Remove the % character from the beginning of lines 250-266 and add the % 
character to the beginning of lines 270-274 and lines 278-280. MATLAB 
ignores lines of code beginning with the % character during function execution. 

To implement Runge-Kutta method of order two 

Remove the % character from the beginning of lines 270-274 and add the % 
character to the beginning of lines 250-266 and lines 278-280. 

To implement Euler ’s method 

Remove the % character from the beginning of lines 278-280 and add the % 
character to the beginning of lines 250-266 and lines 270-274. 

Starting the A VDS /MATLAB Simulation 

The following provides an example of running the F-16 piloted simulation. A similar 
procedure is used for the other provided simulations. 

To start the piloted F-16 simulation 

Step 1. Start AVDS 
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Step 2. Load the simulation initialization file F16.sim.ini 

a. Select File from the top menu 

b. Select Simulation Init, then Open 

c. Choose the file F16.sim.ini located in the F-16 directory of the simulation 
software package 

Step 3. Load the 3-D aircraft image F-16 (Lo-Res) or F-16 (Hi-Res) 

a. Select Initialize from the top menu 

b. Select Aircraft Image 

c. Select either F-16 (Lo-Res) or F-16 (Hi-Res) 

Step 4. In AVDS, switch to simulation mode and start the simulation 

a. Press SIMULATION button on main toolbar 

b. Press START button on main toolbar 
Step 5. Run AVDS_Matlab_F16 in MATLAB 

a. Set the current directory to the F-16 directory of the simulation software package 

b. Type AVDS_Matlab_F16 on the MATLAB command line followed by ENTER 
To end the simulation 

Press the STOP button located on main toolbar in AVDS 
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V. SPECIFIC AIRCRAFT SIMULATION DETAILS 

Details specific to each aircraft simulation package are given below. Unique 
simulation characteristics such as aerodynamic and engine models, aircraft controls, and 
valid ranges for the states and controls are given. Angle of attack, sideslip angle, and 
control surface deflection inputs to the aerodynamic models are in degrees, unless 
specified otherwise. Aircraft angular velocities are in radians per second, unless 
specified otherwise. For a list of the files included in each package, refer to tables A-l 
through A-8 in the Appendix. 


F-16 SIMULATION 

The General Dynamics F-16 is a single-seat, multi-role fighter with a blended 
wing / body and a cropped delta wing planfonn with leading edge sweep of 40°. The 
wing is fitted with leading edge flaps and trailing edge flaperons (flaps / ailerons). Tail 
surfaces are swept and cantilevered. The horizontal stabilator is composed of two all- 
moving tail plane halves, while the vertical tail is fitted with a trailing edge rudder. 
Thrust is provided by one General Electric F110-GE-100 or Pratt & Whitney F100-PW- 
220 afterburning turbofan engine mounted in the rear fuselage. 

The aircraft was modeled with controls for 8th, 8 e , S a , and 8 r . Aerodynamic force and 
moment data were derived from low-speed static and dynamic forced oscillation wind 
tunnel tests conducted on a 16% scale model of the F-16 flown out of ground effect, with 
landing gear retracted, and no external stores 3 ’ 9 . Static aerodynamic data are in tabular 
form as a function of angle of attack and sideslip over the ranges- 10° < a< 45° and 
-30° < j3< 30°. Dynamic data is provided in tabular fonn at zero sideslip angle over the 
angle of attack range -10° < a < 45°. Each non-dimensional aerodynamic force and 
moment coefficient is built up from component functions. Aerodynamic coefficients are 
referenced to a center of gravity location at 0.35 c . Corrections to the flight center of 
gravity position are made in the coefficient build-up equations. 

The aerodynamic database was simplified slightly by dropping second order 
dependencies (e.g., dependence of longitudinal aerodynamic forces on sideslip angle). 
Throttle deflection is limited to the range 0 < <$/, < 1 , elevator deflection is limited to 
-25° < 8 e < 25°, aileron deflection is limited to -21.5° < 8 a < 21.5°, and rudder deflection 
is limited to -30° < 8, < 30°. 

Engine thrust data is in tabular form as a function of power level, altitude, and Mach 
number over the ranges 0 ft < h < 50,000 ft and 0 < M < 1, for idle, military, and 
maximum power settings. Engine power level dynamics and gyroscopic effects 9 are 
included for this aircraft simulation, in the manner described previously. Figure 7 shows 
a three-view of the F-16. 


29 





Figure 7. Three-View of the F-16 


Aerodynamic Model and Database 

The non-dimensional aerodynamic force and moment coefficients for the F-16 vary 
nonlinearly with flow angles (a, /?), aircraft angular velocities (p, q, r), and control 
surface deflections (S e , S a , S r ). Moment coefficients C n and C m include a correction for 
the center of gravity position. The coefficients are computed as follows: 
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Data tables for the component functions of the non-dimensional coefficients are 
defined in the file F16_aero_setup.m. The non-dimensional aerodynamic force and 
moment coefficients are computed for this simulation in F16_aero.m, using linear 
interpolation in the functions cxo.m, czo.m, cmo.m, clo.m, cno.m, dlda.m, dldr.m, dnda.m, 
dndr.m, and dampder.m. 


Engine Model and Database 

The F-16 is powered by a single afterburning turbofan jet engine, which is modeled 
taking into account throttle gearing and engine power level lag. The engine angular 
momentum h eng is assumed to act along the aircraft X body-axis with a fixed value of 
160 slug-ft 2 /s. Throttle gearing is implemented in the file tgear.m, which outputs 
commanded engine power in percent of full power [0,100] using throttle setting input 
[0,1]. Commanded power P c is a function of 5,h, given by: 


J64.944 if S th <0.77 

[217.38^-117.38 if 8 th >0.77 


Engine power level dynamic response is modeled using a first order lag with the time 
constant computed by the function rtau.m. The module pdot.m calculates the rate of 


31 



change of power level with time using the difference between P c and P a . The actual 
power level derivative P a is given by: 


K=^-(Pc-p„) 

^ eng 

where 
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60 
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Engine thrust is modeled as a function of h, M, and P a . Data tables for idle, military, 
and maximum power as a function of h and M are defined in the file 
F16_engine_setup.m. Linear interpolation and thrust calculation is perfonned in the 
module F16_engine.m. Thrust is computed from: 
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Mass Properties 

The mass properties of the F-16 are given in the file F16_massprop.m. These 
properties include the weight of the aircraft, x c . g ., lx, h, Iz, and Ixz ■ This module also 
calculates the constants c 1 ,c 2 ,...,c 9 for use in the body-axis moment equations. The 

longitudinal center of gravity location x c . g . is given in fraction of mean aerodynamic 
chord, and can be adjusted within this module to account for different longitudinal center 
of gravity locations. Table 1 lists the mass properties used in the F-16 simulation. 


Table 1. Mass Properties of the Simulated F-16 


Parameter 

Weight 

(lb) 

lx 

(slug-ft 2 ) 

h 

(slug-ft 2 ) 

(slug-ft 2 ) 

Ixz 

(slug-ft 2 ) 

Value 

20,500 

9,496 

55,814 

63,100 

982 


F-106B SIMULATION 

The Convair F-106B is a two-seat all-weather interceptor incorporating a delta wing 
planform with a single swept vertical tail surface. Elevon control surfaces, which 
perform both elevator and aileron control functions, are located at the trailing edges of 
the wing. The rudder control surface is located at the trailing edge of the vertical tail. 
Thrust is provided by a single Pratt & Whitney J75-P-17 afterburning jet engine mounted 
in the rear fuselage. 

The aircraft simulation has controls for S t h, 8 e , 8 a , and 8,. Low-speed static and 
dynamic forced oscillation wind tunnel tests on a 15% scale model of the aircraft 7 
provided aerodynamic force and moment data. Tests were conducted in the NASA 
Langley 30- by 60-foot wind tunnel, with the model flown out of ground effect, with 
landing gear retracted, and no external stores 10 ’ 11 . Static and dynamic aerodynamic data 
are in tabular fonn as a function of angle of attack over the range -5° < a< 50°. Each 
non-dimensional aerodynamic force and moment coefficient is built up from component 
functions. Aerodynamic coefficients are referenced to a center of gravity location at 
0.275 c , and corrected to the flight center of gravity position in the coefficient build-up 
equations. 
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Throttle deflection is limited to the range 0 < c),/, <1 , elevator deflection is limited to 
-28° < S e < 28°, aileron deflection is limited to -7° < 8 a < 7°, and rudder deflection is 
limited to -25° < 8, < 25°. 

The F-106B engine has maximum thrust values of 24,500 lb and 17,200 lb with and 
without afterburner, respectively. A scaled version of the F-16 engine model was used 
for this simulation. Idle, military, and maximum power settings were scaled using the 
maximum thrust for the F-106B engine. Throttle gearing, gyroscopic effects, and engine 
power level dynamics were carried over from the F-16 engine model. Figure 8 shows a 
three-view of the F-106B aircraft. 



JL JL 



Figure 8. Three-View of the F-106B 


Aerodynamic Model and Database 


The non-dimensional aerodynamic force and moment coefficients for the F-106B 
vary nonlinearly with flow angles (a, ft), aircraft angular velocities (p, q, r ), and control 
surface deflections (8 e , 8 a , 8). Moment coefficients C n and C m include a correction for 
the center of gravity position. The coefficients are computed as follows: 
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Data tables for the component functions of the non-dimensional coefficients are 
defined in the file FI 06B _aero _setup.m. The non-dimensional aerodynamic force and 
moment coefficients are computed for this simulation in F106B_aero.m, using linear 
interpolation in the functions cdo.m, clo.m, cmo.m, dldb.m, dlda.m, dldr.m, dndb.m, 
dnda.m, dndr.m, dydb.m, dydr.m, and dampder.m. 

Engine Model and Database 

Thrust scaling for the F-106B engine is applied in the function 
FI 06B _engine_setup.m. Linear interpolation and thrust calculation is performed in the 
module F106B_engine.m. Functions tgear.m, pdot.m, and rtau.m are unaltered from the 
F-16 engine simulation. 

Mass Properties 

Aircraft mass properties are given in the module F 106B jnassprop.m. Table 2 lists 
the mass properties used in the F-106B simulation. 
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Table 2. Mass Properties of the Simulated F-106B 


Parameter 

Weight 

(lb) 

lx 

(slug-ft 2 ) 

h 

(slug-ft 2 ) 

Iz 

(slug-ft 2 ) 

Ixz 

(slug-ft 2 ) 

Value 

29,776 

18,634 

177,858 

191,236 

5,539 


F-14 SIMULATION 

The Grumman F-14 is a two-seat, carrier-based, multi-role fighter with a variable 
geometry wing and leading edge sweep varying between 20° and 68°. Small canard 
surfaces extend from the leading edges of the fixed portion of the wing to provide extra 
longitudinal stability. Spoilers are located on the upper surfaces of the wing and provide 
aileron control during low speed flight. Trailing edge flaps are located at the rear portion 
of the wing and extend over most of the span, while slats are fitted to the wing leading 
edges. Two all-moving horizontal tail surfaces are located at the rear of the fuselage and 
provide aileron and elevator control using differential and symmetric deflections, 
respectively. Twin vertical tails, also located at the rear of the fuselage, are fitted with 
trailing edge rudder control surfaces. Speed brakes deflect symmetrically from a location 
at the rear fuselage between the bases of the vertical tails. Two Pratt & Whitney TF30-P- 
412A afterburning turbofan engines, located at the rear of the fuselage, provide thrust for 
the aircraft. 

The aircraft was modeled in the unswept wing configuration, with controls for S,/„ 8 e , 
8 a , 8, 8 sp , and 8 S /,. Aircraft aerodynamic data was taken with the vehicle flown out of 
ground effect, with landing gear retracted, and no external stores 12 ’ 13 . Static 
aerodynamic data are in tabular form as a function of angle of attack, sideslip angle, and 
Mach number over the ranges -5° < a < 50°, -20° < j3 < 20°, and 0 <M< 1. The 
aerodynamic effects of angular rates are also in tabular form over the same angle of 
attack range. Each non-dimensional aerodynamic force and moment coefficient is built 
up from component functions. Aerodynamic coefficients are referenced to a center of 
gravity location at 0.16 c , and corrected to the flight center of gravity position in the 
coefficient build-up equations. 

Throttle deflection is limited to the range 0 < <$* < 1 , elevator deflection is limited to 
-30° < 8 e < 12°, aileron deflection is limited to -21° < 8 a < 21°, rudder deflection is 
limited to -30° < 8, < 30°, wing spoiler deflection is limited to 0° < 8 sp < 55°, and speed 
brake deflection is limited to 0° < S s b ^ 60°. 

Tabular data for idle, military, and maximum thrust were used to model the thrust 
produced by the F-14 engine. Engine dynamics and throttle gearing were carried over 
from the F-16 engine model; however, engine angular momentum was neglected. 
Figure 9 shows a three-view of the F-14. 
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Aerodynamic Model and Database 

The non-dimensional aerodynamic force and moment coefficients for the F-14 model 
vary nonlinearly with relative airflow quantities ( a . ; (A M), aircraft angular velocities (p, 
q, r ), and control surface deflections (S e , S a , S r , S sp , S s b). Moment coefficients C n and C m 
include a correction for the center of gravity position. The coefficients are computed as 
follows: 


Cx=C Xo (a,/3) + AC Xs + AC X 


sb 



( sp 

+ AC X 

cj 

1 

■60° 

<•<< 

’V 55 ° 

f 55 J 


C Y = C y (a,p)+C Y ( ct.p)8 a + C r ( a,fi)S r + C Y ( a)5 sp 

u °a °r °sp 

C Y (a) p + C Y (cc )r 


+ 


f b ' 


\ 2v tj 

C z =C z (a,p) + AC z „ + AC Z 


°sb 




+ AC z 

fV 

+ 

r qc ' 

■60° 

V 60y 

\ P =^ 

V55y 


wJ 


Cz q («) 


37 



sp 


C,=C h (a.p) + c, (a,p)S„ +C, S _ (a,p)S r + C, s _ (a,p)r lwJ (M)S l 
\C, (a) P + C, (a) r 


sp 


+ 


f b ' 


v2 V t j 


C m =C m (a,/3) + AC + AC m 


S sb=™° 


V 60 y 


+ AC„ 


’V 55 ° 


'V 

V55y 


V sp ,m{M) + [^-]cm (a) 
\ lv t) 1 


+ ( X, 


c.g.,ref x c.g. ) ^ Z 


)Cz 


c„ = C„„ (a,p)+c ns (a./S) 5„ + C ns ( a,p)S r + C (a)S sp + 

u sp 


f b_ ' 
v2 V tJ 


C n p ( a ) P + C n> . (a)r 


be ' 

I I ^ n rr vnC 


c.g.,ref A c.g. ) ^ Y 


)Cy 


where 


AC X , S =C Xs (a,jB,S e )-C Xo (a,jB) 


AC 


x, 


C 


S sb=™° 




' S S b= 60 ' 


(a)-C x (a) 


AC X =C X (a)-C x (a) 

■S =55° X, 8 =55° ° V 7 

sp sp 


AC z, s = c z, Sg (a, fi,S e )~ Cz 0 {a. P) 


AC 


z. 


z., ,J a )- c z n (a) 

sb 


= C 

S sb = 60° ^ =60° ' ' 


zlC 


z. 


= C 


z,, 


/^- c z n («) 


V 55 ° ’ S sp =55 ° ~° 

=C m -(a.p.S e )-C ma {a.p) 


% ~ m ’8 e 


AC„ 


AC, 


C, 

C, 


m,„ J a )-C mi ,(a) 

} sb~ 


"S sb = 60° ""<U=60° ' 1 0 


w „ w /«, , „ ^m 0 (bx) 

S sp = 55 S sp = 55 


rj sp i (M) = wing spoiler rolling moment efficiency factor 
given as a function of Mach number 

Vsp, m ( M ) = wing spoiler pitching moment efficiency factor 
given as a function of Mach number 
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Data tables for the component functions of the non-dimensional coefficients are 
defined in the file F14_aero_setup.m. The non-dimensional aerodynamic force and 
moment coefficients are computed for this simulation in F14_aero.m, using linear 
interpolation in the functions etamsp.m, etalsp.m, dndr.m, dnda.m, dlda.m, dldr.m, 
dczsp.m, dczsb.m, dczds.m, dcxsp.m, dcxsb.m, dcxds.m, dcmsp.m, dcmsb.m, dcmds.m, 
czo.m, cyo.m, cydsp.m, cydr.m, cyda.m, cxo.m, cno.m, cndsp.m, cmo.m, clo.m, cldsp.m, 
and dampder.m. 

Engine Model and Database 

Data tables for the F-14 engine are defined in the function F 14 _engine_setup.m. 
Linear interpolation and thrust calculation is perfonned in the module F14 engine. m. 
Functions tgear.m, pdot.m, and rtau.m are unaltered from the F-16 engine simulation. 

Mass Properties 


Mass properties for the aircraft are given in the module F14_massprop.m. Table 3 
lists the mass properties used in the F-14 simulation. 

Table 3. Mass Properties of the Simulated F-14 


Parameter 

Weight 

(lb) 

lx 

(slug-ft 2 ) 

h 

(slug-ft 2 ) 

Iz 

(slug-ft 2 ) 

I X7. 

(slug-ft 2 ) 

Value 

48,669 

58,500 

227,000 

276,000 

-2,820 


F-4 SIMULATION 

The McDonnell Douglas F-4 is a two-seat all-weather fighter / bomber. The aircraft 
is fitted with a low-mounted swept wing with wingtip dihedral. The tail section consists 
of an all-moving horizontal stabilator placed in a cathedral configuration and a single 
vertical tail. The trailing edge of the main wing houses control surfaces for the aileron 
and flaps. The trailing edge of the vertical tail accommodates a rudder control surface. 
Thrust is provided by two General Electric GE-J19-17 afterburning jet engines mounted 
at the rear of the fuselage. 

The simulation has controls for 8th, 8 e , S a , and 8,. Aerodynamic data were acquired 
with the model flown out of ground effect, with landing gear retracted, and no external 
stores 14 . Static and dynamic aerodynamics are modeled in terms of polynomial 
functions, which limit control surface effects to first order dependencies, angle of attack 
static dependencies to second order, and sideslip angle static dependencies to third order. 
Polynomial fits for each non-dimensional aerodynamic force and moment coefficient are 
valid over an angle of attack range of -15° < a < 55° . Aerodynamic coefficients are 
referenced to a center of gravity location at 0.289 c , and corrected to the flight center of 
gravity position in the coefficient build-up equations. 
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Throttle deflection is limited to the range 0 < c),/, <1 , elevator deflection is limited to 
-21° < 8 e < 7°, aileron deflection is limited to -15.5° < 8 a < 15.5°, and rudder deflection 
is limited to -30° < 8, < 30°. A scaled version of the F-16 engine model is used with this 
simulation. Thrust for idle, military, and maximum power settings were scaled according 
to the maximum thrust values of the F-4 power plant. Throttle gearing and engine power 
level dynamics were carried over from the F-16 engine model, neglecting engine angular 
momentum effects. Figure 10 shows a three-view of the F-4 aircraft. 



Figure 10. Three-View of the F-4 


Aerodynamic Model and Database 

The non-dimensional aerodynamic force and moment coefficients for the F-4 model 
vary nonlinearly with flow angles (a, /3), aircraft angular velocities {p, q, r ), and control 
surface deflections (8 e , 8 a , 8,). Moment coefficients C„ and C m include a correction for 
the center of gravity position. The coefficients are computed as follows: 
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For a< 15° 


C x =-0.0434 + 2.39xl0“ 3 a + 2.53xl0“V 2 -107xl0“ 6 ct^ 2 +9.5xl0“ 4 4 
-8.5xl0“ 7 £/? 2 +f^^l(8.73xl0“ 3 +0.001^-1. 75xl0“ 4 er) 

l *2 V t J* > 

C Y =-0.012/3 + \.55x\0~ 3 5 r -%x\0~ 6 8 r a 

^ 1 QAL A 

+ (2.25xl0“ 3 ^ + 0.0117r-3.67xl0“ 4 m +1.75xl0“ 4 r£ e ) 

K 7l2V t J 

C z =-0.131-0.0538a-4.76xl0“ 3 4-3.3xl0“ 5 4a-7.5xl0“ 5 ^, 2 

+ f^l(_ 0 . 111 + 5 . 17 x 10 -3a- 1 . 1 x 10 -3« 2 ) 

q = -5.9Sx\0~ 4 /3 - 2.83x\0~ 4 aj3 + L5lx\0~ 5 a 2 j3 

-8 a (6.1xl0“ 4 +2.5xl0“ 5 a-2.6xl0“ 6 a 2 ) 

-S r (-2.3xl0“ 4 +4.5xl0“ 6 a) 

^ 1 QA/ 

+ — - (- 4.l2xl0~ 3 p- 5.24xl0~ 4 pa + 4.36xl0~ 5 pa 2 
y 7T 2V t J 

+ 4.36xl0“ 4 r + 1.05xl0“ 4 m + 5.24xl0“ 5 r4) 

C m =-6.61xl0“ 3 -2.67xl0“ 3 a-6.48xl0“ 5 ^ 2 
-2.65xl0“ 6 a/? 2 -6.54 x 10“ 3 4 -8.49xl0“ 5 4a 
+3.74xl0“ 6 ^ 2 -3.5xl0“ 5 ^ 2 

+ f^rl ( _ 0.0473 - 1.57x1 0 -3 dr) + [x c g re f - x cg ) C z 
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C n =2.28x10“ 3 ^ + 1.79x10”V 3 +1-4x10“ 5 ^ 

+7.0x1 0“ 6 8 cl a- 9.0x1 0“ 4 ^,+4.0x1 0“ 6 8 r a 

^ 1 QAA A 

+ (-6.63xl0“ 5 /?-1.92xl0“ 5 /;a +5.06xl0“V« 2 

K 7l2V t J 

- 6. 06 x 1 0“ 3 r - 8. 73 x 1 0“ 5 rS e + 8. 7 x 1 0“ 6 rd e a) 

~ ~b] ^ Xc ' 8 ' ,re f ~ Xc ' 8 ' ^ 

For 15° < a< 30° 

C x = 0.141-0.0154a + 2.96xl0“V-3.72xl0“ 4 /? 2 
+ 4.14xl0~ 5 a/? 2 - 9.12 x10“V/? 2 +1.82x10“ 3 4 

-7.3xl0“ 5 £,a + -0.0602 + 2.04xl0“ 3 a) 

l *2v t / 

C Y =-2.08xl0“ 2 ^ + 6.07xl0“ 4 a^ + 2.37xl0“V 3 
-3.64xl0“ 7 a^ 3 +2.3xlO“ 3 ^-5.9xlO“ 5 ^a 

^1 QA A 

+ (- 1.62x1 0 -3 /? + 3. 32x1 0 _4 /?a 

V n 2 v t) 

+ 0.031 Ir-1.4xl0“ 3 ra + 1.75xl0“ 4 r4) 

C z =-0.608- 0.022a - 6.77xl0“ 3 4 + 9.7xl0“ 5 4a-7.5xl0“ 5 4 2 

+ f l 80 ^] (i.i36 -0.1418a + 3.11xl0“ 3 a 2 ) 

l ^2 V t J 
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q = -1.29xl0“ 2 ^ + 1.04xl0“ 3 a^-2.02xl0“ 5 a 2 ^ 

+ 1.36xl0“ 5 ^ 3 -1. 13xl0“ 6 a^ 3 +2.01xl0“ 8 a 2 ^ 3 
-4(7.74xl0“ 4 -1.9xl0“ 5 a)- ( 5, (-2.0xl0“ 4 +5.0xl0“ 6 a) 

^ -t QAL 

+ (2.78xl0“ 3 /7-2.79xl0“ 4 /7«-6.81xl0“ 3 r 

[x2V t J 

+ 6. 46 x 1 0“ 4 ra + 5 . 24 x 1 0“ 5 ) 

C m = 0.0549-6.08xl0“ 3 a-1.69xl0“ 4 /? 2 + 5.64xl0“ 7 a/? 2 
- 8.14xl0“ 3 4+l.lxl0“ 4 4a-3.5xl0“ 5 4 2 

+ (-0-0951 + 1.4x10 ~ 3 a) + (x c ^ ref -x c . g )C z 

C n =1.02xl0“ 2 ^- 5. 12xl0“ 4 a^- 5.27x10“ V 3 

+ 3.79xl0“ 7 «A 3 + 9.1xl0“ 5 ^+3.0xl0“ 6 <5 a « 
-1.37xl0“ 3 ^.+3.8xl0“ 5 ^.a 

+ (0.0236/? -2.5xl0“ 3 /?a + 6.25xl0“ 5 /?a 2 

K z 2 v t) 

+ 6.2xl0“ 4 r-4.89xl0“ 4 ra-8.73xl0“ 5 r4 + 8.7xl0“ 6 r4a) 



( x c.g.,ref x c.g. 
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For cc > 30° 

C x =-0.0326 -2. 16xl0“ 3 a + 4.89xlO“V - 1.24xl0' 4 /? 2 
+ 1.076 x10“ 5 q^ 2 -1.54x10“V^ 2 +7.5x10“ 4 4 

-3.9x 10“ 5 <W* + (-0.026 + 8.73xl0“ 4 a) 

l ^ V t J 

C Y =-2.095xl0~ 3 jS-6.36xW~ 5 ajS-2.15xl0~ 5 jS 3 
+5.42xl0~ 7 ajS 3 +1.4xlO~ 3 S r -2.6x10 ~ 5 S,.a 

+ (0. 196/7 - 9.27 xlO -3 /7a + 1.0 lx 10 _4 /?a 2 

U 2 V t J 

+0. 0326r - 2. 55 x 1 0“ 3 ra + 3. 26 x 1 0“ 5 ra 2 + 1 . 75 x 1 0 -4 rS e ) 
C z =-0.891-0.01 146a + 6.2xl0“ 3 4 - 5.4xl0“ 3 £ e a 
+ 6.2xl0“ 6 4a 2 -7.5 x10“ 5 4 2 

+ (0.589- 0.0494a + 6.11xlO“V) 

l J 

Q =1.18xl0~ 2 jS-5.29xl0~ 4 ajS + 4.88xl0~ 6 a 2 jS 
-2.2x1 0 -5 jS 3 +9.05x10 ~ 7 ajS 3 -9.08x10 ~ 9 a 2 jS 3 
-5.0xl0~ s S a -S r (-9.0xl0“ 5 +1.8 x10“ 6 «; 

^ -t QA/a^ 

+ (-0.0428/7 + 1. 82xl0“ 3 /7a-1.94xl0“ 5 pa 2 

U 2 V t ) 

+0.073r-3.02xl0“ 3 ra + 3.14xl0“ 5 ra 2 +5.24xl0“ 5 /^ e ) 

C m =7.3xlO“ 3 -5.5xlO“ 3 a-7.93xlO“ 3 (5 e + 8.23xlO“ 5 <5 e a 

-3.5x 10“ 5 4 2 +fl^9£_l (0. 16 -0.0 101a + 1.05xl0 _4 a 2 ) 
U 2V t J 

+ { x c.g.,ref ~ x c.g. ) Q 
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C„ =-9.23xl0“V + l-52xl0“ 4 a y 5 + 1.62xl0“V 3 
-3.46xl0" 7 o^ 3 + 1.5xl0^^ a -6.8x10^^ 
+1.2xl0“ 5 ^«-<5 e (1.67xl0“ 4 A-5.56xl0“ 6 «A 


-3.81xl0“ 7 ^ 3 +1.27xl0“ 8 a^ 3 ) 


+ 


^ 1807? A 


(0.0385/? - 1.73x1 0 _3 /?a + 1.92xl0 _5 /?a 2 


\*W tJ 

-0. 0202r + 3. 67 x 1 0“ 4 ra + 2. 5 8 x 1 0“ 3 rS. 


( c ' 


- 1 . 3 1 x 1 0“ 4 rd e a + 1 . 69 x 1 0“ 6 rS e a 2 ) - [-j (x cgref - x cg ) C y 


The non-dimensional aerodynamic force and moment coefficients are computed for 
this simulation in F4_aero.m. 

Engine Model and Database 

Scaled thrust values for the F-4 power plant are applied in the function 
F4 _engine_setup.m. Linear interpolation and thrust calculation are performed in the 
module F4_engine.m. Functions tgear.m, pdot.m, and rtau.m are unaltered from the F-16 
engine simulation. 

Mass Properties 


Aircraft mass properties are given in the module F4_massprop.m. Table 4 lists the 
mass properties used in the F-4 simulation. 

Table 4. Mass Properties of the Simulated F-4 


Parameter 

Weight 

(lb) 

lx 

(slug-ft 2 ) 

h 

(slug-ft 2 ) 

Iz 

(slug-ft 2 ) 

I xz 

(slug-ft 2 ) 

Value 

38,924 

24,970 

122,190 

139,800 

1,175 


FASER SIMULATION 

A radio-controlled aircraft, called the NASA Langley FASER, has the form of a 
high-wing general aviation aircraft. The aircraft was modeled with controls for 5th, 5 e , 5 a , 
and 5 r . Non-dimensional static and dynamic aerodynamic data were estimated 15 using 
airfoil data and approximations based on aircraft geometry, for the aircraft out of ground 
effect. The aerodynamic data is in tabular form as a function of angle of attack, Mach 
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number, lift coefficient, and thrust over the ranges -10° < a < 30°, 0.05 <M< 0.3, 
- 1 < C L <1.6, and 0 lb < T < 10 lb. Each non-dimensional aerodynamic force and 
moment coefficient is built up from component functions, which are provided in the 
stability-axis frame. Aerodynamic coefficients are referenced to a center of gravity 
location at 0.273 c , and corrected to the flight center of gravity position in the coefficient 
build-up equations. 

Throttle deflection is limited to the range 0 < 4 a < 1 , elevator deflection is limited to 
-25° < 4 < 25°, aileron deflection is limited to -25° < 4 < 25°, and rudder deflection is 
limited to -30° < S, < 30°. FASER is fitted with a single electric motor that drives a 
propeller to generate thrust. Figure 1 1 shows a three-view of the FASER. 




Figure 11. Three-View of the FASER 


Aerodynamic Model and Database 

The non-dimensional aerodynamic force and moment coefficients for the FASER 
vary nonlinearly with relative airflow quantities {a, j3, M), lift coefficient, {Cl ), thrust 
(7), aircraft angular velocities {p , q, r ), and control surface deflections (4, 4, 4). 
Moment coefficients C n and C m include a correction for the center of gravity position. 
The coefficients are computed as follows: 


C x = C L sin 


( 0 , 71 ^ 

vl 80 y 


-C D cos 


/ 

vl 80 y 
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where 


Cy=C r/j (M) 


( jBn^ 
vl80y 


+ Cy s (M) 


( 8 ,.n^ 


V 180 y 


+ 


v2 V tJ 


Cy (cr ,M^) p s + Cy (a,M)i 


C 7 = —Cj cos 


Ps 


f an ^ 


V 180y 


■C D sin 


f an ^ 


V 180y 


C l = C l,s cos 


f an ^ 


V 180y 


' C n,S sin 


f an ^ 


V 180y 


^ qc^ 


C m =C (a, M)+ i— c (a, M) + C ( a,M ) 

yzv t ) 1 °e 

+ AC m + [ x c.g.,ref — ^c.g. ) 


V 180 y 


Q = c /i5 sin 


an 

l80j 


+ Q,S my ^ j( X c.g.,n?/ x c.g.)^Y 


C L =C L (aM) + Cr ( a,M ) 



+ AC, + 

f qc ^ 

1 180 y 

L, ,j~’ 

wJ 


C Lq (a,M) 


C D =C DQ {a,M) + AC Dq 


C ls =C lp {M,C L ) 


f fin ^ 

vl80y 




+ V M WS (a ’ M) 




V 180 y 


+ 


V 2 ^y 


Q (M)p s +Cj (a,M,C L ) 


V Ps 


r s 


( j 8 n^ 


C n , s =C np {a,M)\^yC n/; (M,C l ) 


V 180 y 


+ C n§ (a,M) 


( 8 ,-n\ 


V 180 ) 


+ 


r b_' 

yW tJ 


c ' (a,M,Cj^+C„ (a,M,C L ) 


Ps 


r s 


AC l , t =C L ' T (a,M ,T)-C Lq (a,M) 
AC D ' T =C DT (a,M,T)-C Do {a,M) 
AC m<T = C Mt {a,M.T)-C {a.M) 
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Ps =PCOS 



+ r sin 



r s 


-p sin 



+ r cos 



Data tables for the component functions of the non-dimensional coefficients are 
defined in the file FASER_aero_setup.m. The non-dimensional aerodynamic force and 
moment coefficients are computed for this simulation in FASER_aero.m, using linear 
interpolation in one-, two-, and three-dimensional interpolation functions intl.m, int2.m, 
and int3.m. 

Engine Model and Database 

Propulsion for the FASER is modeled neglecting engine angular momentum. 
Throttle gearing is implemented in the file tgear.m, which outputs commanded engine 
power in percent of full power [0,100] using throttle setting input [0,1]. Commanded 
power P c is a linear function of S t i„ given by: 

p c (s, h )=ms lh 

Engine power level dynamic response is modeled using a first order lag with time 
constant computed by the function rtau.m. The module pdot.m calculates the rate of 
change of power level with time using the difference between P c and P a . The actual 
power level derivative P a is given by: 

Pa=—(P C -Pa) 

^ eng 

where 

1.0 if (P c -P„)< 25.0 

0.2 if >50 0 

10 -f^][(^-/;,)- 2 5.0] if 25.0 < {P c - P a ) < 50.0 

Engine thrust is modeled as a linear function of P a . Thrust calculation is perfonned 
in the module FASER_engine.m using the equation: 

T = 0AP a 
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Mass Properties 

Mass properties of the aircraft are given in the module FASERjnassprop.m. Table 5 
lists the mass properties used in the FASER simulation. 


Table 5. Mass Properties of the Simulated FASER 


Parameter 

Weight 

(lb) 

lx 

(slug-ft 2 ) 

h 

(slug-ft 2 ) 

Iz 

(slug-ft 2 ) 

I xz 

(slug-ft 2 ) 

Value 

16.5 

1.324 

1.946 

0.627 

0.317 


HL-20 SIMULATION 

The NASA HL-20 is a manned lifting-body spacecraft. It is designed to be placed 
into earth orbit by an expendable launch vehicle and re-enter the atmosphere for an 
unpowered approach and landing. Aircraft control surfaces consist of an all moving 
vertical tail (rudder), two upper and lower body flaps, and two trailing edge wing flaps. 
Elevator control is provided by symmetric deflection of the wing flaps, aileron control by 
differential deflection of the wing flaps, positive flap control by symmetric deflection of 
the lower body flaps, negative flap control by symmetric deflection of the upper body 
flaps, and differential flap control by differential deflection of the upper and lower body 
flaps. 

The HL-20 was modeled with controls for 8 e , S a , S r , S pl , 8 n f, Sdf, and 8/ g . Low-speed 
(M <0.6) wind tunnel data for a scale model of the aircraft flown out of ground effect 
with landing gear retracted and no external stores were taken in the Langley 30- by 60- 
foot tunnel and Calspan 8-foot transonic tunnel 16 ’ 17 . Landing gear effects were based on 
scaled Space Shuttle data, and dynamic derivative data were found using estimation 
techniques. Static and dynamic aerodynamic data are modeled with polynomial fits as a 

function of angle of attack and sideslip angle over the ranges - 1 0° <a< 30° and 
-10° < J3< 10°. Each non-dimensional aerodynamic force and moment coefficient is built 
up from component functions. Aerodynamic coefficients are referenced to a center of 
gravity location at 0.54 c , and corrected to the flight center of gravity position in the 
coefficient build-up equations. 

Elevator deflection is limited to -30° < S e < 30°, aileron deflection is limited to 
-30° < 8 a < 30° , rudder deflection is limited to -30° < 8, < 30°, positive flap deflection 

is limited to 0° < 8 p f < 60°, negative flap deflection is limited to -60° < 8 n f < 0°, 
differential flap deflection is limited to -60° < 8 d f < 60°, and landing gear deflection is 
limited to 0° < Si g < 98°. An engine model is not included in this simulation since the 
HL-20 uses no propulsion during atmospheric flight after re-entry. Figure 12 shows a 
three-view of the HL-20. 
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s 


Figure 12. Three-View of the HL-20 


Aerodynamic Model and Database 

The non-dimensional aerodynamic force and moment coefficients for the HL-20 
model vary nonlinearly with flow angles (a, j3), aircraft angular velocities (p , q, r), and 
control deflections (8 e , 8 a , 8 r , S p f, S nh 8 d f, Si g ). Moment coefficients C n and C m include a 
correction for the center of gravity position. The coefficients are computed as follows: 


c* = C X(i + C Xg S e + C Xs % I +c Xf I s r I + c Xs s pf + c Xg S„ f + C x |5, 

e a r pj nj uj 


Cy = -0.01242/? + C 7j 8 a + 8 r + C Yg 8 df 


C 7 =C 7 +C 7 S e + C 7c 8 nf + C /c . 8 nf + AC 7 s 
z z 0 z S e e z 8 pf PJ z 5 nf n ? z ’% 
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Q (J ( a - £> ) + (■„ s 3 a + C„ s £>■ + ^ng $df + 




where 

C x = -(7 .362 x 1 0“ 2 - 2 .56 x 1 0 ~ 4 a -2.208x1 0 ~ 4 a 2 - 2.262 x 1 0 “V + 2 .966 x 1 0“ V 
-3.64xl0‘V+9.388xl0“ 1 V’-5.299xl0“ 4 |/?|-4.709xl0“V 2 +8.572xl0“ 5 fi 3 
-4. 199x1 0 -6 /? 4 + 1 .295 x 1 0~ 4 a \f }\ ) 

C x§ = -( - 1 .854 x 1 0“ 4 + 2 .83 x 1 0 “ 6 a -6.966x1 0“ V +1.323x1 0“ V -2.758x1 0~ V ) 

C x = -(9.776 x 1 0 -4 - 2.703 x 1 0“ 5 a - 8.303 x 1 0 ~ V + 6.645 x 1 0“ V - 1 .273 x 1 0“ V ) 

°a 

C x = -(5.812xl0“ 4 + 1.41xl0“ 5 tf-2.585xl0“ 6 er +3.051xl0“V-8.161xl0“V) 

C Y = 3.357xl0“ 3 -1.661xl0“ 5 a-3.28xl0“ 6 er +5.526xlO“V-3.269xlO“ 10 a 4 
C Y = 1.855xlO“ 3 + 1.128xl0“ 5 tf + 6.069xl0“ 6 a 2 -1.78xlO“V -1.886xl0“ 12 a 4 

C F =2.672xl0“ 5 -3.849xl0“ 5 a + 4.564xl0“ 7 er +1.798xlO“V -4.099xl0“ 10 a 4 

5 df 

C 7 =-(- 9.025 xl0“ 2 + 4.07xl0“ 2 a + 3.094xl0“ 5 er +1.564xl0“ 5 a 3 -1.386xl0“ 6 a 4 
z o 

+2.545 xlO -8 a 5 -1.189xlO“ 10 etr 6 + 2.564xl0“ 3 |/?| + 8.501xl0“ 4 /? 2 
-1.156xl0“ 4 /? 3 +3.416xlO“V 4 -4.862xl0“ 4 a|^|; 

C z§ =-(5.140xl0“ 3 +3.683xl0 _5 a-6.092xl0 -6 a 2 +2.818xlO“V -2.459xl0“ 9 a 4 ) 
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c x = -(1.31xl0“ 4 + 1.565xl0“ 6 a 2 -1.542xlO“V) 

5 pf 

C Xx = -(-4.415xl0“ 4 -4.056xl0“ 6 a-4.657xl0“V) 

V 

C x e =-(-6.043xl0“ 4 -1.858xl0“ 5 a + 8xl0“V-4.848xl0“V+1.36xl0“V) 

s df 

AC x,s lg = c x,s lg (<*, S lg ) - C Xq ( a ) 

C z = -(3.7779xl0“ 3 -7.017xl0“ ? a 2 + 1.4xlO“ 10 a 4 ) 

5 pf 

C z = -(3.71 lx 10 -3 - 3 .547 x 1 0 -5 a - 2 .706 x 1 0 -6 a 2 + 2 .93 8 x 1 0 -7 a 3 

V 

-5.552x10 “ 9 a 4 ) 

Ac z ,s lg =c z , Sig ( a ,s lg )-c Zo ( a ) 

C, = 2.538xlO“ 3 +1.963xlO“ 5 a-3.725xlO“ 6 a 2 +3.539xlO“ 8 a 3 -1.778xlO“ 10 a 4 

S a 

C, = 2.26 x 1 0“ 4 - 1 .299 x 1 0 ~ 5 a + 5.565x1 0“ 6 a 2 - 3 .3 82 x 1 0“ V + 6.46 1 x 1 0“ V 

L 0„ 


C. =- 7.453xl0“ 4 -1.811xl0“ 5 a-1.264xl0“V+9.972xl0“V-2.684xl0“V 
V 


C, = -0.526 1357 + 0.029749438a - 0.006960559a 2 + 0.00053564a 3 - 1 .35043x 10“ 5 a 4 

1 


C, = 0.497793859-0. 028742272a + 0. 016480243a 2 -0.001439044a 3 +3.56179xl0 _5 a 4 

V 
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C m = 2.632 xlO -2 - 2.226 xlO“ 3 a-l. 8594 xlO“ 5 a 2 + 6.001xlCrV + 1.828xl0“ 7 a 4 
m o 

-9.733x 10“V + 1.71x 10“ 10 a 6 - 5.233x 10“ 4 \j3\ + 6.795 x 10 ~ 5 jB 2 
-1.993xl0“ 5 ft +1.341xl0“V 4 +6.061xl0“ 5 a|/?| 

C = -1.903xl0“ 3 -1.593xl0“ 5 a + 2.611xl0“ 6 a 2 +5.116xlO“V -1.626xlO“V 

m 8 e 

c ; =-9.896xl0“ 4 -1.494xl0“ 9 a 2 +6.303xl0“ n a 4 

S Pf 

C ; =-1.086xl0“ 3 +1.57xl0“ 5 a + 4.174xl0“ 7 a 2 -1.133xlO“ 7 a 3 +2.723xlO“V 

S nf 

AC m,8 lg = C m,S lg (aAg ) - C m Q («) 

C m = -0.195552417- 0.006776550a + . 00285193a 2 -0.000184146a 3 

q 

+2.45653xl0“ 6 a 4 

C lo = -2.769 xlO -3 - 4.377 x 10 _5 a + 9.952x 10 _6 a 2 -3.642xlO“V +4.692xlO“V 

S a 

C„ c =-1.278xl0“ 3 +1.32xl0“ 5 a-4.72xl0“ 6 a 2 +2.371xl0“ 7 a 3 -3.340xl0“ 9 a 4 

S r 

C„ =-5.107xl0“ 5 +1.108xl0“ 5 a-1.547xl0“ 8 a 2 -1.552xl0“ 8 a 3 +1.413xlO“ 10 a 4 

s df 

C n = 0.385 147809 - 0.00 135606 1 a - 0.00209489a 2 +8.52445 x 10“ 5 a 3 + 1 .02309x 10“ 7 a 4 

p 

C n = -0. 790855552-0. 010770449a-0. 00095452a 2 +0. 000151569a 3 -5.13541xl0 _6 a 4 

"r 

Data tables required for the component functions of the non-dimensional coefficients 
are defined in the file HL20 aero _setup.m. The non-dimensional aerodynamic force and 
moment coefficients are computed for this simulation in HL20_aero.m. Linear 
interpolation is done in the functions cno.m, cmlgear.m, cxlgear.m, and czlgear.m. 

Dynamic Equations 

Thrust effects and the power level state derivative were removed from the dynamic 
equations for the HL-20 simulation. 

Mass Properties 

Aircraft mass properties are given in the module HL20 jnassprop.m. Table 6 lists the 
mass properties used in the HL-20 simulation. 
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Table 6. Mass Properties of the Simulated HL-20 


Parameter 

Weight 

(lb) 

lx 

(slug-ft 2 ) 

h 

(slug-ft 2 ) 

Iz 

(slug-ft 2 ) 

Ixz 

(slug-ft 2 ) 

Value 

19,100 

7,512 

33,594 

35,644 

0 


X-31 SIMULATION 

The NASA / DARPA X-3 1 is a thrust-vectored, single-seat research aircraft designed 
to explore highly maneuverable flight. The aircraft has a cranked delta-wing planform, 
all-moving canard surfaces mounted near the nose, and a single vertical tail. The wing 
houses leading and trailing edge flaps. Trailing edge flaps act as elevon control surfaces, 
which function as separate elevator and aileron control surfaces within the simulation. 
Speed brakes are attached to the rear section of the fuselage and deflect symmetrically 
outward when actuated. A single General Electric F404-GE-400 non-afterburning 
turbofan engine provides thrust for the aircraft. Three paddles are symmetrically 
attached at the rear exit of the engine nozzle to provide thrust-vectoring capabilities. 
However, thrust vectoring is not implemented in this simulation. The effect of symmetric 
outboard deflection of the three thrust-vectoring paddles is modeled as an additional 
speed-braking device. 

The aircraft was modeled with controls for 5th, 4 , 5 a , 5 r , 5 C , S s b, and 4 » Low-speed 
static wind tunnel tests provided force and moment aerodynamic data for a 19% scale 
model of the aircraft 18 . Tests were conducted in the Langley 14- by 22-foot subsonic 
tunnel with the model flown out of ground effect, with landing gear retracted, and no 
external stores. Dynamic data were obtained from forced oscillation wind tunnel 
testing 19 . Static aerodynamic data are in tabular fonn as a function of angle of attack and 
sideslip angle over the ranges -5° < a< 55° and -18° < j3< 18°. Each non-dimensional 
aerodynamic force and moment coefficient is built up from component functions. 
Aerodynamic coefficients are referenced to a center of gravity location at 0.5286 c , and 
corrected to the flight center of gravity position in the coefficient build-up equations. 

Throttle deflection is limited to the range 0 < 5 t h < 1, elevator deflection is limited to 
-40° < 5 e < 30°, aileron deflection is limited to -30° < 5 U < 30°, rudder deflection is 
limited to -30° < 5, < 30°, canard deflection is limited to -20° < 5 C < 60°, speed brake 
deflection is limited to 0° < 5 s h < 50°, and thrust vectoring paddle deflection is limited to 
0° < 5 tv < 50°. Data from the F-16 engine model was scaled using the maximum thrust 
values for the X-31 engine. Throttle gearing and engine power level dynamics were 
carried over from the F-16 engine model; however, engine angular momentum was 
ignored. Figure 13 shows a three-view of the X-31. 
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Figure 13. Three- View of the X-31 


Aerodynamic Model and Database 

The non-dimensional aerodynamic force and moment coefficients for the X-3 1 model 
vary nonlinearly with flow angles (a, ft), aircraft angular velocities (p, q, r), and control 
surface deflections (8 e , 8 a , 8 r , 8 C , 8 s b, 8 h ). Moment coefficients C„ and C m include a 
correction for center of gravity position. The coefficients are computed as follows: 
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Data tables for the component functions of the non-dimensional coefficients are 
defined in the file X3 1 _aero _setup.m. The non-dimensional aerodynamic force and 
moment coefficients are computed for this simulation in X31_aero.m, using linear 
interpolation in the functions ca_can.m, dca_dele.m, dca_dsb.m, dcadtv.m, 
cnorm_can.m, denorm dele. m, dcnormdsb.m, dcnormdtv.m, cy_beta.m, dcy_drdr.m, 
dcy_dail.m, dcy_dsb.m, dcy_dtv.m, cm_can.m, dcm_dele.m, dcm_dsb.m, dcm_dtv.m, 
cl_beta.ni, dcl drdr.m, dcl dail.m, dcl_dsb.m, dcl_dtv.m, cn_beta.m, dcn drdr.m, 
dcn_dail.m, dcn_dsb.m, dcn_dtv.m, caq.m, cyr.m, cyp.m, cnormq.m, clr.m, clp.m, cmq.m, 
cnr.m, and enp.m. 

Engine Model and Database 

Thrust scaling for the X-31 engine is applied in the function X3 1_engine_setitp.ni. 
Linear interpolation and thrust calculation is perfonned in the module X31 engine. m. 
Functions tgear.m, pdot.m, and rtau.m are unaltered from the F-16 engine simulation. 

Mass Properties 

Aircraft mass properties are given in the module X31 jnasspr op. m. Table 7 lists the 
mass properties used in the X-3 1 simulation. 
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Table 7. Mass Properties of the Simulated X-31 


Parameter 

Weight 

(lb) 

lx 

(slug-ft 2 ) 

h 

(slug-ft 2 ) 

Iz 

(slug-ft 2 ) 

Ixz 

(slug-ft 2 ) 

Value 

16,000 

3,553 

50,645 

49,367 

-156 


A- 7 SIMULATION 

The Vought A-7 Corsair II is a single-seat tactical attack fighter based on the F-8 
Crusader platform. The aircraft is controlled with ailerons located on the trailing edge of 
the main wing, a rudder control surface located at the trailing edge of the vertical tail, and 
an all-moving horizontal stabilator located at the rear of the fuselage. Thrust is provided 
by a single non-afterburning Pratt & Whitney TF30-P-408 turbofan engine mounted in 
the rear fuselage. 

Aircraft controls are provided for 5th, 4, 5 a , and 5 r . Low-speed static wind tunnel 
tests and dynamic forced oscillation wind tunnel tests provided force and moment data 
for a scale model of the aircraft flown out of ground effect, with landing gear retracted, 
and no external stores 20 . Static and dynamic aerodynamic data are in tabular form as a 
function of angle of attack and sideslip over the ranges 0° < a< 90° and 0° < j3< 90°. 
Each non-dimensional aerodynamic force and moment coefficient is built up from 
component functions. Aerodynamic coefficients are referenced to a center of gravity 
location at 0.35 c , and corrected to the flight center of gravity position in the coefficient 
build-up equations. 

Throttle deflection is limited to the range 0 < 4, <1 , elevator deflection is limited to 
-25° < 5 e < 15°, aileron deflection is limited to -25° < 5 a < 25°, and rudder deflection is 
limited to -30° < 5, < 30°. The A-7 engine has a maximum thrust of 13,400 lb. A scaled 
version of the F-16 engine model was used for this simulation. Idle and military power 
thrust data were scaled according to the maximum thrust for the A-7 engine. Throttle 
gearing, engine gyroscopic effects, and engine power level dynamics were carried over 
from the F-16 engine model. Figure 14 shows a three-view of the A-7 aircraft. 
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Figure 14. Three-View of the A-7 


Aerodynamic Model and Database 

The non-dimensional aerodynamic force and moment coefficients for the A-7 vary 
nonlinearly with flow angles (a, /?), aircraft angular velocities {p, q, r), and control 
surface deflections (S e , S a , S r ). Moment coefficients C n and C,„ include a correction for 
the center of gravity position. The coefficients are computed as follows: 
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Data tables for the component functions of the non-dimensional coefficients are 
defined in the file A7_aero_setup.m. The non-dimensional aerodynamic force and 
moment coefficients are computed for this simulation in A7_aero.m, using linear 
interpolation in one- and two-dimensional interpolation functions intl.m and int2.m. 


Engine Model and Database 

Thrust scaling for the A-7 engine is applied in the function A7 _engine_setup.m using 
the F-16 engine data, omitting afterburner data. Throttle gearing, implemented in 
tgear.m, outputs commanded engine power in percent of full power [0,100] using throttle 
setting input [0,1]. Commanded power P c is a linear function of 5th, given by 


P,(S„) = 100<5„, 


Functions pdot.m and rtau.m are unaltered from the F-16 engine simulation. 

Linear interpolation and thrust calculation is performed in the module A7_engine.m, 
scaling F-16 engine data to the maximum thrust of the A-7 engine. Scaled thrust is 
computed from 


T = 


" 13,400 " 

v 12,680 y 


(r a „+o.oi(r„„ 



Mass Properties 

Aircraft mass properties are given in the module A 7 jnasspr op. m. Table 8 lists the 
mass properties used in the A-7 simulation. 
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Table 8. Mass Properties of the Simulated A-7 


Parameter 

Weight 

(lb) 

lx 

(slug-ft 2 ) 

h 

(slug-ft 2 ) 

Iz 

(slug-ft 2 ) 

I xz 

(slug-ft 2 ) 

Value 

22,699 

16,970 

65,430 

76,130 

4,030 
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VI. GENERATING 3-D MODELS FOR AVDS 

Three-dimensional aircraft models used within AVDS are constructed using a 
separate CAD drafting program. Models created with CAD are converted to an input 
format suitable for AVDS. AVDS supplies a utility that allows conversion of Drawing 
Interchange Format (DXF) 3-D aircraft models into AVDS format. The following 
sections provide a tutorial for creating and adding a 3-D aircraft model into AVDS for 
use in AVDS / MATLAB piloted simulations. 


MODEL CONSTRUCTION REQUIREMENTS 

Three-dimensional aircraft models must be constructed in AutoCAD 8 DXF format 
using surface modeling techniques. The AVDS conversion program does not recognize 
models created using solid modeling techniques. Surface models have exteriors 
composed of meshes made of three and four sided polygons whose vertices are located in 
three-dimensional space. The total number of polygons in a surface mesh detennines 
model complexity. AVDS / MATLAB interactive simulations are computationally 
intensive; therefore, it is necessary to create 3-D models that are as simple as possible in 
order to minimize graphics rendering requirements. Complex models with dense surface 
meshes can slow AVDS / MATLAB interactive simulations considerably on less 
powerful computers. 

Many 3-D CAD programs support both surface modeling and the DXF fonnat. Ruled 
surfaces, tabulated surfaces, and edge surfaces are used extensively during surface model 
creation. A surface created using any of these tools is recognized by the AVDS 
conversion utility. However, revolved surfaces and two-dimensional lines are not 
recognized during conversion. 

Aircraft 3-D models should be created using two-dimensional surfaces for the wings 
and horizontal / vertical tails. This reduces the complexity of the aircraft surface mesh 
while not excessively detracting from aircraft realism. The fuselage, canopy, and other 
significant aircraft structures should be created as simply as possible using low-density 
surface meshes. Following these guidelines will simplify model making as well as 
prevent any unnecessary slow-downs during piloted simulations run on less powerful 
computers. 


MODELING EXAMPLE OF A SIMPLE 3-D AIRCRAFT 

An example is given to illustrate the surface model creation of a simple 3-D aircraft, 
shown in figure 15. It is assumed that the reader has familiarity with using and drawing 
in AutoCAD Release 14 (or higher) or an equivalent 3-D CAD program. For more 
information on CAD design or surface modeling consult any AutoCAD Release 14 (or 
higher) reference book and surface-modeling tutorial such as Reference [21]. Aileron, 
rudder, and elevator control surfaces will be included in the model, while wheels and any 
external stores are ignored. 
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Figure 15. Three-View of the MiG-15 Aircraft 

It is convenient to model aircraft using available three-views such as the one shown 
in figure 15. These drawings may be found on various internet sites, and in references 
such as [22]. Dimensions used during initial modeling of the aircraft should be those of 
the drawing rather than the actual airplane. This simplifies model-making by eliminating 
the need to scale the dimension of every point on the model. Once modeling is complete, 
the finished 3-D aircraft must be scaled to the dimensions of the actual aircraft. 
Dimensions must be in feet for proper implementation into AVDS. The scaling feature is 
available in most CAD programs. The scaled model must also be aligned along the 
world coordinate system (WCS), prior to conversion into the AVDS format, with the 
nose of the model pointing along the positive X-axis and positive Z-axis pointing out the 
top of the model. The origin of the WCS must be positioned at the location of the aircraft 
center of gravity. AVDS uses this origin location as placement for the body, earth, and 
stability axes on the model during interactive simulations. 

The fuselage for this particular aircraft is a relatively simple shape and can be 
modeled as a cylinder with the front and back openings used as the intake and exhaust. 
Figure 16 shows the preliminary front and rear openings of the aircraft. 
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Figure 16. Front and Rear Fuselage Openings 


The AutoCAD command RULESURF is used to create a surface mesh between the 
intake and exhaust. The completed fuselage is shown in figure 17. 
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Figure 17. Modeled MiG-15 Fuselage Section 


The wing, tail, and control surfaces of the aircraft are modeled using planar surfaces. 
This is accomplished using the AutoCAD command 3DFACE. Figure 18 shows the 
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aircraft with the various surfaces added. RULESURF was used to create surface meshes 
for the modeled intake and exhaust. 



The canopy can be modeled as an elongated dome. A wire frame is constructed 
outlining the general shape of the canopy, shown in figure 19. Surface meshes are drawn 
over the wire frame using the RULESURF command. 
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Figure 19. Model with Wire Frame Canopy 
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The finished canopy is shown in Figure 20. The final drawing must be scaled (in 
feet) to the dimensions of the actual aircraft and aligned accordingly with the WCS 
specified previously. 



Figure 20. Completed 3-D MiG-15 Model 


IMPORTING A COMPLETED 3-D CAD MODEL INTO AVDS 

A completed model saved in DXF format is converted to the AVDS format with the 
DXF to AVDS conversion utility. For usage instructions, refer to Appendix E of the 
AVDS User’s Manual . The box marked Two sided polygons should be checked in the 
interface window of the conversion utility prior to conversion, in order to create surface 
meshes composed of front and back sided polygons. If this option is not chosen, meshes 
composed of either front or back sided polygons will be shown in AVDS. This would 
cause surfaces to disappear unexpectedly when viewing the aircraft from certain angles. 

The DXF to AVDS conversion utility has difficulty implementing entire aircraft 
models into the AVDS format. Therefore, it is necessary to separate each major model 
element (e.g., elements that will share the same color, control surfaces elements, canopy, 
etc.) into individually saved DXF files. These individual files are converted separately 
and assembled into a single AVDS craft file by the user. The user must identify section 
color, wing and tail ribbon location, pilot view, control surfaces, and control surface 
rotation axes for each newly compiled AVDS craft file, following the format given in 
Appendix A of the AVDS User’s Manual 2 . 

Control surfaces for new aircraft models must also be identified within the files 
ACname.sim.ini and AVDS_Matlab_ACname.m, where ACname is the aircraft name (e.g. 
F16.sim.ini, AVDS_Matlab_F16, etc.), using the same control surface category 
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designations as were used in the AVDS craft file. Refer to Section 4.2 in the AVDS 
User’s Manual for instructions on creating / modifying ACname.sim.ini files. Figure 21 
shows the input screen for files of this type. Within this configuration window, FCS 
must be set to Basic FCS, AIRCRAFT should be set to MATLAB AC, and Simulation 
Input Type should be set to Positions. Elements within the MATLAB vectors 
InputVector and OutputVector, which have been declared by the user in 
AVDS_Matlab_ACname.m, must also be declared to AVDS via the ACname.sim.ini file. 
The user must declare InputVector elements under the MATLAB AC sub-heading 
section of INPUTS, while OutputVector elements must be declared under the 
Environment sub-heading section under INPUTS. Control surfaces are identified in 
AVDS as sfcOl to sfc!5 under the Environment sub-heading section under INPUTS. 



Figure 21. Input Screen for ACname.sim.ini 
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VII. CONCLUDING REMARKS 

Software has been developed in MATLAB for six-degree-of-freedom nonlinear 
simulations of various aircraft. Trim, linearization, numerical integration, batch 
simulation, 3-D visualization, and real-time piloted simulation software are included for 
all simulations. 

Aerodynamic models were developed using data gathered from open literature reports 
documenting wind tunnel tests, flight tests, and analytic estimation techniques. Most of 
the engine models contained in the aircraft simulation packages were based on the F-16 
engine model, using scaled values for thrust. Other models were used for aircraft engines 
when suitable data was available or when the jet engine model did not apply (e.g., 
FASER and the HL-20 aircraft). 

Interactive piloted simulations can be run using the nonlinear aircraft simulation 
packages and visual flight simulator software. The nonlinear simulations described in 
this report can be used for education, control law development, dynamic analysis, and 
many other purposes. 

The structure of the simulations was designed for easy use and to allow modifications 
with relatively little effort. The collection of aircraft simulations gives the user a range of 
aircraft dynamics to explore. This compilation is also a testament of the applicability of 
the basic simulation structure to a variety of aircraft. 

Future improvements to this work, outside of adding more aircraft, could include 
such changes as adding landing gear dynamics, allowing for more general aircraft e.g. 
locations, including a more realistic gravitational model that accounts for earth geodesy, 
including a non-stationary atmosphere (e.g., turbulence, micro-bursts, winds, etc.), and 
adding aeroelastic effects. Since AVDS currently allows only four control inputs, an 
improvement that allows the user to command a larger number of controls would be 
desirable for interactive simulation. It is anticipated that network capability will be built 
into future editions of the simulations, allowing multiple users to fly interactively with 
(or against) one another in real time. Finally, an obvious improvement for nearly all the 
aircraft simulations would be the addition of a flight control system. This was 
intentionally omitted in the current work, to encourage the use of these aircraft nonlinear 
simulations for control system design and analysis. 
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IX. APPENDIX 

Table A-l. F-16 Simulation Package (Directory name: F-16) 


F-16 Directory 

ab3.m 

F16_deq.m 

atm.m 

F16_engine.m 

AVDS_Matlab_F16.m 

F16_engine_setup.m 

clo.m 

F16_hires.txt 

cmo.m 

F16_lores.txt 

cno.m 

F16_massprop.m 

cnvrg.m 

F16_trm.m 

cxo.m 

gen_F16_model.m 

czo.m 

grad.m 

dampder.m 

icftrm.m 

dlda.m 

Inze.m 

dldr.m 

mksqw.m 

dnda.m 

pdot.m 

dndr.m 

rk2.m 

F16.m 

rk4.m 

F16.sim.ini 

rtau.m 

F16_aero.m 

solve. m 

F16_aero_setup.m 

tgear.m 
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Table A-2. F-106B Simulation Package (Directory name: F-106B) 


F-106B Directory 

ab3.m 

F106B aero.m 

atm.m 

F106B_aero_setup.m 

AVDS Matlab F106B.m 

F106B_deq.m 

cdo.m 

F106B_engine.m 

clo.m 

F106B_engine_setup.m 

cmo.m 

F 1 06B_massprop.m 

cnvrg.m 

F106B_trm.m 

dampder.m 

gen_F106B_model.m 

dlda.m 

grad.m 

dldb.m 

ic_ftrm.m 

dldr.m 

Inze.m 

dnda.m 

mksqw.m 

dndb.m 

pdot.m 

dndr.m 

rk2.m 

dydb.m 

rk4.m 

dydr.m 

rtau.m 

F106B.m 

solve. m 

F106B.sim.ini 

tgear.m 

F106B.txt 
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Table A-3. F-14 Simulation Package (Directory name: F-14) 


F-14 Directory 

ab3.m 

dcmsp.m 

F14_deq.m 

atm.m 

dcxds.m 

F14_engine.m 

AVDS Matlab F14.m 

dcxsb.m 

F14_engine_setup.m 

cldsp.m 

dcxsp.m 

F14_massprop.m 

clo.m 

dczds.m 

F14_trm.m 

cmo.m 

dczsb.m 

gen_F14_model.m 

cndsp.m 

dczsp.m 

grad.m 

cno.m 

dlda.m 

ic_ftrm.m 

cnvrg.m 

dldr.m 

Inze.m 

cxo.m 

dnda.m 

mksqw.m 

cyda.m 

dndr.m 

pdot.m 

cydr.m 

etalsp.m 

rk2.m 

cydsp.m 

etamsp.m 

rk4.m 

cyo.m 

F14.m 

rtau.m 

czo.m 

F14.sim.ini 

solve. m 

dampder.m 

F14.txt 

tgear.m 

dcmds.m 

F14_aero.m 


dcmsb.m 

F 1 4_aero_setup.m 



Table A-4. F-4 Simulation Package (Directory name: F-4) 


F-4 Directory 

ab3.m 

F4_trm.M 

atm.m 

gen_F4_model.m 

AVDS_Matlab_F4.m 

grad.m 

cnvrg.m 

icftrm.m 

F4.m 

Inze.m 

F4.sim.ini 

mksqw.m 

F4.txt 

pdot.m 

F4_aero.m 

rk2.m 

F4_deq.m 

rk4.m 

F4_engine.m 

rtau.m 

F4_engine_setup.m 

solve. m 

F4_massprop.m 

tgear.m 
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Table A-5. FASER Simulation Package (Directory name: FASER) 


FASER 1 

Directory 

ab3.m 

grad.m 

atm.m 

ic_ftrm.m 

AVDS Matlab FASER. m 

inti .m 

cnvrg.m 

int2.m 

FASER. m 

int3.m 

FASER.sim.ini 

lat_aero_tables.mat 

FASERaero.m 

Inze.m 

FASER_aero_setup.m 

lon_aero_tables.mat 

FASERdeq.m 

mksqw.m 

FASER_engine.m 

pdot.m 

FASER_hires.txt 

rk2.m 

FASER_lores.txt 

rk4.m 

FASER massprop.m 

rtau.m 

FASER_trm.m 

solve. m 

genFASER model, m 

tgear.m 


Table A-6. HL-20 Simulation Package (Directory name: HL-20) 


HL-20 Directory 

ab3.m 

HL20.txt 

atm.m 

HL20_aero.m 

AVDS Matlab HL20.m 

FIL20_aero_setup.m 

cmlgear.m 

HL20_deq.m 

cno.m 

HL20_massprop.m 

cnvrg.m 

FIL20_trm.m 

cxlgear.m 

icftrm.m 

czlgear.m 

Inze.m 

gen_HL20_model.m 

mksqw.m 

grad.m 

rk2.m 

HL20.m 

rk4.m 

HL20.sim.ini 

solve. m 


73 





Table A-7. X-31 Simulation Package (Directory name: X-31) 


X-31 Directory 

ab3.m 

dcadtv.m 

icftrm.m 

atm.m 

dcl_dail.m 

Inze.m 

AVDS Matlab X31.m 

dcldrdr.m 

mksqw.m 

caq.m 

dcl_dsb.m 

pdot.m 

cacan.m 

dcldtv.m 

rk2.m 

clp.m 

dcm_dele.m 

rk4.m 

clr.m 

dcmdsb.m 

rtau.m 

cl_beta.m 

dcm_dtv.m 

solve. m 

cmq.m 

dcnormdele.m 

tgear.m 

cm_can.m 

dcnorm_dsb.m 

X31.m 

cnormq.m 

dcnormdtv.m 

X31 .sim.ini 

cnorm_can.m 

dcn_dail.m 

X31.txt 

cnp.m 

dcndrdr.m 

X31_aero.m 

cnr.m 

dcn_dsb.m 

X31_aero_setup.m 

cnvrg.m 

dcndtv.m 

X31 deq.m 

cn_beta.m 

dcy_dail.m 

X31_engine.m 

cyp.m 

dcydrdr.m 

X31_engine_setup.m 

cyr.m 

dcy_dsb.m 

X31_massprop.m 

cybeta.m 

dcydtv.m 

X31_trm.m 

dca_dele.m 

gen_X31_model.m 


dca_dsb.m 

grad.m 



Table A-8. A-7 Simulation Package (Directory name: A-7) 


A-7 Directory 

ab3.m 

gen_A7_model.m 

atm.m 

grad.m 

AVDS_Matlab_A7.m 

icftrm.m 

A7.txt 

intl.m 

A7.m 

int2.m 

A7.sim.ini 

Inze.m 

A7_aero.m 

mksqw.m 

A7_aero_setup.m 

pdot.m 

A7_deq.m 

rk2.m 

A7_engine.m 

rk4.m 

A7_en g i n esetu p . m 

rtau.m 

A7_massprop.m 

solve. m 

A7_trm.m 

tgear.m 

cnvrg.m 
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